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SECTION  I 

INTRODUCTION  AND  SUMMARY 

1 . 1 BACKGROUND 

The  process  by  which  projectiles  penetrating  fuel  cells 
transfer  kinetic  and  chemical  (explosive)  energy  to  the  fuel 
cell  structure  is  termed  hydrodynamic  ram.  Hydrodynamic  ram 
has  proven  to  be  a major  combat-related  threat  to  modern  air- 
craft. The  development  of  design  principles  for  fuel  cells 
less  vulnerable  to  ram  damage  is  a priority  task.  Hence,  a 
number  of  research  progreims  have  been  conducted  to  investigate 
the  coupling  of  projectile  energy,  energy  to  penetrated 
fluids , the  coupling  of  fluid  momentum  to  containing 
structures , and  the  response  of  fuel  cell  structures  to 
fluid  loading^^"^'  7-9)^ 

The  present  program  is  addressed  to  hydrodynamic  ram 
induced  by  high-velocity  fragment  simulators.  When  such  a 
projectile/fragment  penetrates  a filled  fuel  cell  structure, 
energy  is  transferred  to  the  structure  through  pressure  in  the 
fluid.  The  pressure  in  the  fluid  is  generated  by  two  different 
mechanisms,  the  impact  induced  shock  waves,  and  the  s'lbsonic 
flow  ciround  the  penetrating  projectile.  Each  mechanism  will  be 
reviewed  in  detail. 


1.1.1 


Shock  Induced  Pressure 


The  portion  of  the  projectile  trajectory  during  which 
shock  waves  are  generated  is  sometimes  referred  to  as  the  shock 
or  entrance  phase  of  hydrodynamic  ram.  When  the  projectile 
strikes  the  entrance  panel,  a shock  wave  is  driven  into  the 
panel  material.  At  the  panel/fluid  boundary,  this  shock  is 
transmitted  into  the  fluid.  If  the  projectile  is  supersonic 
in  the  fluid  after  panel  perforation,  shock  waves  may  continue 
to  propagate  away  from  the  projectile.  The  pressure  propagates 
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into  the  fluxd  as  an  expanding  guasihemispherical  shock  wave. 
The  amplitude  of  the  wave  decays  rapidly  due  to  geometric 
spreading  and  the  generation  of  release  waves  at  the  entrance 
site.  As  the  shock  wave  propagates  along  the  entrance  panels 
the  panel  is  exposed  to  extremely  high  pressures  for  a short 
time.  The  resulting  panel  motion  can  produce  localized  stress 
states  near  the  entrance  hole  which  are  sufficient  to  initiate 
tensile  cracks  and  gross  tearing  of  the  entrance  panel . 

The  maximum  possible  amplitude  of  the  impact  induced 
shock  in  the  fluid  is  that  which  would  result  from  a one 
dimensional  impact.  The  magnitude  of  the  one  dimensional 
pressure  can  be  calculated  from  the  Hugoniot  relationship 

P = P u U * (1) 

o 

cind 

p U = p,  (U-u)  (2) 

0 X 

together  with  the  measured  P-u  Hugoniot s for  the  fluid 
(water)  , projectile^^^^ , and  cell  wall^^^^ , and  the  tradi- 
tional assumptions  that  pressure  and  particle  velocity  are 
continuous  across  interfaces . The  magnitude  of  the  impact 
pressure  computed  from  (1)  and  (2)  is  very  high,  as  illustrated 
in  Figure  1.  For  comparison,  the  tensile  strength  of  aluminum 
is  only  about  500  MN/m^.  The  peak  shock  pressure  in  the  fluid 
will  be  less  than  the  value  shown  in  Figure  1 for  relatively 
pointed  projectiles  (e.g.,  of  impact  radius  less  than  the  wall 
thickness) . Due  to  geometric  speeding,  the  peak  pressure  in 
the  fluid  will  decay  inversely  with  radial  distance  (as  1/r) . 

In  addition,  release  waves  generated  at  the  boundaries  of  the 
impact  region  (entrance  panel)  overtake  the  shock  front  and 
further  reduce  its  intensity. 


*P,  p , p- , u,  and  U denote,  respectively,  the  pressure  across 
the  ^shock  front,  the  density  of  unshocked  material,  the  density 
of  shocked  material,  the  particle  velocity  change  across  the 


shock,  and  the  propagation  velocity  of  the  shock  with  respect 


to  the  undisturbed  material. 
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Figure  1.  Peak  Shock  Pressure  in  Tanks  Prenetrated  by  Blunt 

Steel  Projectiles.  (a)  Walls  of  Negligible  Thick- 
ness. (b)  Thick  Aluminum  Walls. 

1.1.2  Drag  Related  Pressure 

In  the  drag  or  cavity  phase  of  a hydrodynamic  ram 
event,  the  projectile  is  slowed  by  drag  forces  as  it  passes 
through  the  fluid.  The  momentum  imparted  to  the  fluid  by  the 
projectile  produces  a cavity  behind  the  projectile  and  a pres- 


sure pulse  throughout  the  fluid  which  loads  the  cell  panels 
over  a large  area.  If  the  projectile  is  tumbling  or  breaks  up, 
significant  fluctuations  in  pressure  can  result.  Fluid  motion 
and  pressure  usually  persist  even  after  the  projectile  has 
come  to  rest  or  exited  the  tank.  Occasionally,  local  failure 
inside  panels  result  from  the  drag  phase  loading.  More  fre- 
quently, failure  occurs  at  the  fasteners  which  are  unable  to 
resist  the  net  force  on  the  panel. 

A number  of  analytical  treatments  of  this  phase  of 

(7  9) 

hydrodynamic  ram  have  been  developed  ' , The  starting  point 

for  such  analyses  is  the  energy  loss  rate  of  the  projectile 


id 
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given  by 


6x 


(3) 


where  E is  the  kinetic  energy  of  the  projectile,  x is  the  dis- 
tance along  trajectory,  is  the  drag  coefficient,  a is  the 
projectile  cross  sectional  area  and  a is  the  projectile  velocity. 
The  energy  loss  rate  is  assumed  to  be  due  to  the  work  which  the 
drag  phase  pressure  does  on  the  penetrating  projectile. 

1.1.3  Exit  Phase 

The  exit  phase  occurs  when  the  projectile  perforates 

a panel  and  leaves  the  cell.  In  the  exit  phase,  the  exit  panel 

is  first  stressed  by  fluid  which  is  accelerated  and  compressed 

ahead  of  the  projectile  and  then  struck  by  the  projectile  itself. 

The  coupling  of  hydrodynamic  ram  with  ballistic  loading  can 

result  in  extensive  radial  cracking  of  the  exit  panel  in  the 
(3) 

impact  region'  . When  such  an  exit  panel  failure  occurs,  it 
is  generally  attributed  to  a steep  increase  in  fluid  pressure 
as  the  projectile  approaches  the  wall,  followed  immediately  by 
projectile  impact.  In  some  cases  the  high  pressure  preceding 
the  projectile  is  sufficient  to  initiate  radial  cracking  before 
perforation  occurs. 

1 . 2 OBJECTIVES 

The  research  program  described  in  this  report  had  two 
principal  objectives.  The  first  was  to  evaluate  two  available 
analytic  codes  for  prediction  of  hydrodynamic  ram  effects.  A 
second  objective  was  to  obtain  an  understanding  of  ram  pheno- 
menology for  high  velocity  spherical  fragments. 

1.3  APPROACH  OF  PROGRAM 

The  hydrodynamic  ram  prediction  codes  which  appeared  to 
be  most  promising  were  the  BR-lA(HR)  finite  element  code  and 
the  AFTON  finite  difference  code.  (These  codes  are  described 
in  detail  in  Section  V of  this  report) . The  experiments  were 
designed  to  provide  the  input  data  which  the  codes  require  and 
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to  provide  data  with  which  the  output  of  the  codes  could  be 
compared . 

The  principal  target  employed  was  a replica  of  an  A-10 
fuselage  tank  with  a capacity  of  760)1  (200  gal  ) , For  con- 
venience and  safety,  the  tank  was  filled  with  water.  The  panel 
material  was  1.6  :nnm  thick  aluminum.  Two  alloys  were  selected, 
7075-T6  and  2024  T-3.  They  find  widespread  use  in  aircraft  and 
have  well  known  material  properties.  Some  panels  were  backed 
with  AVCO  Thermarest^  ballistic  foam;  this  material  was  selected 
because  previous  work  has  suggested  its  effectiveness  in 

/ g \ 

reducing  ram  damage  , and  because  its  "crush-up''  properties 
have  been  well  determined. 

Early  in  the  progreun  it  was  found  that  important  entrance 
phase  phenomena  occur  at  very  early  times  and  are  highly 
localized.  To  aid  the  study  of  these  effects,  a series  of 
experiments  was  conducted  in  a small  shock  tank  of  28)^  volume. 

It  was  felt  that  the  complications  caused  by  asymmetrical 
entrance  holes  and  tumbling  projectiles  would  defeat  an  attempt 
to  assess  the  basic  validity  of  the  computer  treatments. 
Therefore,  spherical  projectiles  were  used  throughout  this 
program.  The  velocities  employed  spanned  the  range  expected  for 
high  velocity  fragmejits  in  combat,  1.5  to  2.4  km/s.  The  tra- 
jectory was  at  90®  to  the  entrance  panel. 

The  data  required  to  exercise  the  hydrodynamic  ram  codes 
are  the  projectile  position  as  a function  of  time  in  the  fluid, 
and  the  velocity  of  the  shock  front  in  the  fluid.  In  addition, 
pressure  measurements  within  the  fluid  were  obtained.  The  codes 
were  verified  by  comparisons  between  calculated  and  measured 
panel  displacements. 

Measurement  of  panel  out-of-plane  displacement  was 
carried  out  using  a moire  fringe  technique.  This  technique  has 
a temporal  resolution  of  80  us  and  a displacement  resolution 
of  0.3  mm.  Measurements  of  projectile  trajectory  was  determined 
from  a 7000  f/s  framing  camera  record  of  the  progress  of  the 
projectile  through  the  fluid.  Shock  velocity  measurement  was 
made  with  a 10^  f/s  framing  camera.  A number  of  piezoelectric 
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pressure  transducers  were  placed  in  the  tank  to  obtain  pressure 
data. 


1.4  SUMMARY  OF  RESULTS 

Twenty-nine  shots  were  completed  for  which  data  was 
obtained  for  out-of -plane  entrance  or  nide  panel  displacement, 
projectile  drag  coefficient,  cavity  expansion  rate,  shock 
velocity,  and  pressure.  The  NWC  fluid  drag  code'  ' and  a 
fluid  shock  code  were  applied  to  the  experimental  configura- 
tions and  pressure  input  was  generated  for  the  structure  code. 
The  BR-lA(HR)  code  was  adapted  for  the  fluid  loads  on  the 
structure.  The  boundary  conditions  corresponding  to  the 
laboratory  target  tank  were  set  up.  The  California  Research 
and  Technology  (CRT)  AFTON  code  was  also  applied  to  the  experi- 
mental configurations.  The  AFTOIf^NONSAP  calculation  showed  good 
qualitative  and  quantitative  agreement  with  the  experimental 
results.  The  BR-IA(HR)  results  showed  poor  agreement  with 
experimental  measurements.  The  difficulty  was  traced  to  the 
fluid  pressure/wall  motion  coupling  scheme. 

The  experimental  data  support  the  following  statements: 

1.  High  velocity  fragments  caused  severe  damage  to 
entrance  panels. 

2.  In  these  experiments  foam  of  thickness  equal  to 
the  projectile  diameter  effectively  eliminate*.!  entrance  panel 
damage. 

3.  Aluminum  alloy  2024-T3  was  superior  to  7075-T6  for 
defeating  hydrodynamic  ram. 

4.  For  tlie  configuration  and  threat  used  in  these  tests, 
side  panel  failures  were  not  initiated. 

5.  The  mcire  fringe  technique  provided  an  excellent  way 
to  analyze  panel  motion;  distance-time  resolutions  of  0.3  mm 
and  80  us  were  readily  achievable. 

6.  The  NWC  pressure  code  is  a valuable  technique  for 
quantitative  calculation  of  early-time  behavior  and  qualitative 
calculation  of  late-time  events  when  reflections  dominate  the 
pressure. 
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SECTION  II 

EXPERIMENTAL  FACILITY 

Two  ballistic  ranges  at  the  University  of  Dayton  Research 
Institute  were  instrumented  for  study  of  hydrodynamic  ram  effects 
in  fuel  cells.  One  of  those  ranges  was  dedicated  to  the  main 
series  of  experiments  which  consisted  of  measurements  of  ram 
induced  fuel  cell  panel  motion,  fluid  pressure,  and  projectile 
trajectories.  The  second  range  was  devoted  to  study  of  early 
time  shock  effects  in  fuel  cells. 

2.1  MAIN  EXPERIMENTAL  FACILITY 

The  hardware  and  instrumentation  for  the  principal 
range  used  in  this  program  is  shown  in  Figure  2 . 


Figure  2.  Diagram  of  Range  Hardware  and  Instrumentation  for 
Main  Experimental  Facility. 
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It  consisted  of  a projectile  launcher,  velocity  measurement 
stations,  moire  fringe  displacement-measuring  equipment,  and 
an  instrumented  replica  A-10  fuselage  section. 
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2.1.1  Projectiles 

The  projectiles  used  in  this  program  were  spherical 
fragment  simulators  made  from  annealed  ball  bearings  (52100 
chrome  manganese  steel).  The  two  sizes  used  were  5.5  and  11  g, 
with  diameters  of  11.1  and  14.3  mm.  The  microhardness  of  the 
projectiles  was  280  Hv  (Vickers  hardness  scale,  10''  g load)  . 

The  spheres  were  launched  in  serrated  polycarbonate 
sabots  which  were  stripped  by  aerodynamic  forces  and  stopped 
by  metal  baffles  within  the  blast  tank.  Launch  velocities  were 
varied  from  1.4  to  2.3  km/s. 

2.1.2  Launch  Apparatus 

The  projectile  launch  apparatus  consisted  of  a breech, 
barrel,  and  blast  tank.  The  specially-designed  breech  had  a 
chamber  capacity  up  to  110  g.  The  barrel  had  a bore  diameter 
of  20  mm  and  a length  of  2.3  m ; the  muzzle  end  was  inserted 
into  the  blast  tank.  The  blast  tank  was  evacuated  to  a pres- 
sure of  about  10  kN/m‘  (0.1  Atm).  Baffle  plates  within  the 
blast  tank  served  to  contain  the  sabot  as  well  as  to  reduce  the 
air  blast  within  the  range. 

The  velocity  measurement  section  of  the  range  consisted 
of  aluminum  foil-mylar  switches  and  two  flash  x-rays.  The 
switches  were  constructed  of  two  sheets  of  aluminum  foil 
separated  by  a thin  sheet  of  mylar.  The  first  switch  was 
placed  29  cm  downrange  of  the  blast  tank  and  the  second  switch 
was  placed  55  cm  downrange  of  the  first.  In  addition,  for  most 
shots,  a third  switch  was  mounted  on  the  target  tank  front  wall 
1.8  m beyond  the  first  switch.  The  radiographs  were  used 
to  determine  projectile  position  and  integrity  at  the  switches. 
The  first  two  switches,  as  shown  in  Figure  2,  were  connected  to 
passive  pulsers  which  triggered  both  the  flash  x-rays  and  a 
10  MHz  counter- timer . On  several  shots,  a 1 MHz  counter- timer 
was  connected  between  the  second  and  third  switches. 
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To  compute  projectile  velocities,  times  were  ta.cen 
from  the  counter  readings  and  distances  were  computed  in  two 
steps.  First,  the  distance  traversed  by  the  projectile  was 
determined  approximately  using  the  distance  between  the  foil 
switches.  Secondly,  the  distance  was  determined  more  accurately 
from  readings  of  the  projectile  position  in  the  radiogiaphs. 
Projectile  velocity  could  be  determined  to  be  v/ithin  * 0.5% 
using  the  radiographs  and  to  within  * 1.0%  using  only  the 
distances  between  ttie  foil  switches. 

2.1.3  Target  Tank 

The  target  tank  consisted  of  two  sections  of  a replica 
A-10  fuselage.  The  sections  were  permanently  bolted  together 
and  the  center  ports  removed.  The  overall  dimensions  of  the 
tank  are  indicated  in  Figure  3.  A photograph  of  this  tank  is 
shown  in  Figure  4 . 


Figure  3.  Drawing  of  760 e Target  Tank. 
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Fiqure  4.  Photoqraph  of  760C  Taroet  Tank,  with  Entrance  and 
Side  Panels  Removed.  Supports  for  holdinvj  five 
pressure  transducers  are  mounted  independent  of  tank. 

The  tank  was  filled  with  approximately  760C  of  tap  water.  The 
two  bulkhead  panels  and  one  side  panel  were  used  as  test  panels. 
The  test  panels  were  made  of  1 mm  thick  aluminum  alloy  and 
were  fastened  with  9.5  mm  (3/>’.  in)  bolts  at  38  mm  intervals. 

A layer  of  rubber  between  the  panels  and  the  frame  prevented 
leakage.  The  other  side  panel  was  constructed  of  19  mm  thick 
aluminum  with  a 0.3  m by  1.0  m viewing  port.  The  viewing  port 
was  covered  with  a 12  mm  thick  polycarbonate  window.  The  bottom 
of  the  tank  was  made  of  1.6  mm  thick  7075-T6  alumin<im  backed 
with  10  mm  of  plywood.  The  frame  was  9.5  irim  thick  stainless 
steel.  The  top  of  the  tank  was  open.  The  tank  was  orientated 
so  that  the  bulkhead  panels  were  the  entrance  and  exit  panels 
for  the  projectiles. 

2.1.4  Target  Tank  Instrun\entation 

2. 1.4.1  Moird  fringe  apparatus 

A moire  fringe  optical  system  was  mounted  on  a 
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movable,  adjustable  tripod,  and  was  directed  at  cither  the 
entrance  or  the  side  test  panel.  This  system  is  described  in 
detail  in  Section  III. 

The  moire  fringe  data  was  recorded  with  a high- 
speed framing  camera  (Fastax) . The  speed  cf  the  camera  was 
determined  either  by  monitoring  the  drive  sprocket  rotation 
with  a reluctance  loop  and  oscilloscope,  or  by  using  a built-in 
timing  light  source.  Two  speeds  were  used,  6000  and  12,000  t s. 
Illumination  was  accomplished  with  "Press  50"  flash  lamps. 

The  film  employed  was  Eastman  Double  X developed  in  Acufine 
developer  for  10  rain. 

2. 1.4. 2 Projectile  drag  measurements 

The  progress  of  the  projectile  through  the  tank,  and 
the  development  of  the  cavity,  were  recorded  with  a high-speed 
framing  camera.  This  camera  photographed  the  interior  of  the 
tank  through  the  polycarbonate  side  window.  Its  speed  was 
monitored  with  an  oscilloscope,  and  was  either  6000  or  12,000 
f/s.  Illumination  for  this  camera  consisted  of  six  reflector 
photoflood  lamps. 

2. 1.4. 3 Pressure  measurements 

Five  (occasionally  four)  piezoelectric  pressure 
transducers  were  placed  in  the  tank.  They  were  mounted  in 
aluminum  caps  threaded  onto  20  ram  diameter  galvanized  steel 
water  pipe  and  inserted  into  the  tank  from  the  top.  The  tip 
of  the  transducer  protruded  slightly  from  the  cap  to  eliminate 
the  possibility  of  an  air  bubble  forming. 

Most  of  the  transducers  were  PCB  10''A04.  This 
transducer  is  calibrated  to  70  MN/m*  and  has  a rise  time  of 
1 us.  For  measurements  near  the  entrance  point,  a Kistler 
Model  207C3  transducer  was  sometimes  employed.  It  has  a peak 
pressure  rating  of  700  MN/m- . The  transducer  outputs  were 
recorded  on  cathode  ray  oscilloscopes.  The  oscilloscope 
trigger  signal  was  generated  from  the  third  foil  switch 
located  at  the  entrance  panel. 


2.2  SHOCK  TANK  RiU^GE 

The  shock  tank  rango  w^as  assembled  to  obtain  data  tor 
shock  wave  propagation  .ia  the  iifonediate  vicinity  ot  the  entrance 
hole  of  a penetrcited  taci  coll  . A diagram  ot  this  range  is  pre- 
sented in  Figure  5.  The  projectiles,  the  launch  aiJparatus,  and 
the  toil  switches  were  essentially  identical  to  those  used  on 
the  principal  range.  The  target  tank  was  a cube  ot  2 8v'  capacity. 
The  sides  v/cre  6.3mm  thick  plexi<il.uu;  and  tlie  enlraiu-o  panel 
was  1.6  mm  thick  aluminum. 

Instrumentation  consisted  ot  a Beckman  and  Whitley  Model 
300  framing  camera  and  spark  light  source.  The  camera  exposed 
48  frames  at  speeds  up  to  4.5  x 10'’  f/s.  The  camera  runs  con- 
tinually during  the  expcruAcnt  and  data  is  obtained  by  timing 
the  discharge  of  the  spark  light  source.  The  spark  was 
triggered  from  the  third  foil  switch  mounted  near  the  entrance 
panel,  ns  shown  in  Figure  5.  The  camera  framing  rate  was 
generally  6 x 10'  f/s. 


Figure  5.  Ins l ramontation  for  Shock  Tank  Range. 
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SECTION  111 

ANALYSIS  OF  MOIRE  FRINGE  DATA 


Tho  reduction  of  moire  I'rinyc  data  to  yield  inlormatiou 
on  out-of-plane  displacement  constituted  a major  effort  of  this 
procjram.  Two  approaches  were  developed;  empirical  and  theoro- 
ticaJ . In  the  empirical  approach,  the  locations  of  frinyes  in 
space  were  determined  by  means  of  suitable  fiducial  data.  Tho 
position  of  each  frinye  on  the  target  panel  was  then  used  to 
obtain  the  out-of-plane  displacement  of  the  panel  at  that  point. 
In  the  theoretical  approach,  the  form  of  the  fringe  pattern  was 
derived  from  the  geometry  and  dimensions  of  the  optical  system. 
Tho  theoretical  approach  is  more  powerful  and  can  be  automated 
with  computer  programs.  A computer  program  based  on  the 
theoretical  treatment  was  developed  and  debugged  with  the  aid 
of  the  empirical  solutions. 

The  combination  f theoretical  and  empirical  data 
reduction  techniques  provided  a deep  understanding  of  moire 
fringes  and  the  properties  of  the  fringe  patterns  used  in  those 
experiments.  In  the  following  sections  the  moire  apparatus  is 
described  in  detail,  the  theoretical  model  of  fringe  production 
is  explained,  and  the  empirical  data  reduction  technique  is 
described. 


3.1  DESCRIPTION  OF  MOIRE  OPTICAL  APPARATUS 

The  inoir€  optical  apparatus  consisted  primarily  of  two 
Ronchi  rulings,  commercially  available  lenses,  a high  intensity 
light  source,  and  a high  speed  framing  camera.  The  apparatus 
is  shown  schematically  in  Figure  6.  The  optical  elements  were 
mounted  onto  one-meter  long  optical  benches  which  were  bolted  to 
a 1.2  m square,  19  mm  thick  aluminum  table.  The  table  could 
be  both  elevated  and  tilted. 
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The  light  source  chosen  was  a "Press  50"  flash  lamp  with 
an  intensity  of  more  than  0.8  Mlm  for  at  least  30  ms,  as  shown 
in  Figure  7.  This  was  sufficient  light  for  recording  approxi- 
mately 180  frames  of  information  at  a f railing  rate  of  6,000  f/s. 
The  Ronchi  rulings  were  glass  slides  with  evenly-spaced,  finely- 
ruled  lines,  separated  by  a distance  equal  to  the  line  width. 

The  rulingr  had  78.7  lines  per  centimeter  (200  lines  per  inch) . 

A Fas tax  high-speed  framing  camera  was  used  to  record  the 
moire  fringes. 

An  optical  block  diagram  of  the  moire  fringe  apparatus 
is  shown  in  Figure  8.  Light  from  the  flash  lamp  was  collected 
by  the  condensing  lens  and  focused  to  a point  inside  the  pro- 
jection lens.  The  projectile  lens,  a 80  mm-f/2.8  lens,  was 
used  to  project  an  image  of  the  specimen  grid  (a  Ronchi  ruling) 
onto  the  target  wall.  The  image  of  the  specimen  grid  on  the 
tank  wall  was  called  the  target  grid.  Due  to  the  angle  of  pro- 
jection, the  pitch  spacing  between  grid  lines  of  the  target 
was  not  constant  as  illustrated  in  Figure  9. 
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Figure  7.  Intensity  vs.  Time  for  Flash  Light  Source. 
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The  target  grid  was  focused  through  a 75  mm  - f/1.5  lens  onto 
a master  grid  (a  second  Ronchi  ruling) . The  geometric  inter- 
ference of  the  master  grid  with  the  image  of  the  target  grid 
produced  fringes  characteristic  of  the  target  wall  topography. 
The  high  speed  framing  camera  recorded  this  changing  moire 
pattern  as  the  target  was  impacted  by  steel  spheres. 

A simple  equivalent  model,  valid  in  the  region  of  the 
target,  was  used  to  describe  the  moire  optical  apparatus.  In 
the  model,  rays  radiated  from  the  projection  lens  and  converged 
in  the  taking  lens  at  points  on  the  optic  axes  very  near  the 
front  principal  planes  of  the  respective  lenses.  These  points 

.*  A 

were  called  the  moire  source  point  and  the  moire  convergence 
point.  The  straight  lino  connecting  these  points  defined  a 
reference  line  for  the  moire  apparatus  called  the  baseline. 

A triangle  was  defined  by  the  principal  axes  of  the 
projection  and  photographic  optical  elements  and  the  baseline 
between  the  moire  source  and  convergence  points.  This  triangle, 
illustrated  in  Figure  10,  was  called  the  moire  triangle.  To 
facilitate  measurement  of  the  triangle,  the  principal  optical 
axes  were  designed  to  lie  parallel  to  the  center  lines  of  the 
optical  benches. 

3.2  THEORETICAL  ANALYSIS  OF  MOIRE  FRINGES 

Out-of-plane  motion  of  the  target  panel  caused  corres- 
ponding changes  in  the  moire  fringe  pattern.  Measurements  of 
this  pattern  and  the  dimensions  and  orientation  of  the  moire 
triangle  were  used  to  calculate  displacements  of  the  target. 

A geometric  model  of  the  moire  system,  valid  in  the 
region  of  the  wall,  was  constructed  to  facilitate  data  reduc- 
tion. The  model,  illustrated  in  Figure  11,  also  helped  clarify 
the  method  of  fringe  information.  The  model  of  the  optical 
system  was  simplified  by  considering  all  light  to  radiate  from 
a source  point  inside  the  projection  lens.  Similarly,  reflected 
light  from  the  target  converged  to  a point  inside  the  taking 
lens.  Equivalent  specimen  and  master  grids  were  drawn  between 
the  target  and  source  and  the  target  and  convergence  points. 
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Figure  11.  Moire  Geometric  Interference  Diaqraw. 

The  source  and  convergence  points  are  points  at  which  the  rays 
appear  to  eminate  from  or  converge  to.  These  points  are 
physically  located  at  the  lens  principal  plane  furthest  from 
the  grid. 

The  equivalent  specimen  grid  was  positioned  a distance 
€2  in  front  of  the  source  point;  light  and  dark  rays  were 
produced  as  light  was  alternately  passed  and  blocked  by  this 
grid.  This  resulted  in  the  formation  on  the  tank  wall  of  a 
target  grid  with  a varying  pitch.  Light  from  the  target  grid 
was  alternately  passed  and  blocked  by  an  equivalent  master 
grid  located  on  the  photographic  optical  axis,  a distance  f^  in 
front  of  the  convergence  point. 

Regions  in  which  light  was  blocked  by  the  Ronchi  rulings 
are  heavily  shaded  in  the  illustration.  The  heavily  shaded 
areas  of  the  master  grid  are  extended  as  lightly  shaded  areas 
to  facilitate  the  description  of  the  model.  Any  point  of  the 
target  that  falls  in  the  heavily  shaded  region  is  not  visible 


’0 


because  the  light  to  that  point  is  blocked  by  the  specimen 
grid.  Any  point  on  the  target  which  lies  in  a lightly  shaded 
region  is  also  not  visible  because  it  lies  behind  an  opaque 
portion  of  the  master  grid. 

Lines  can  be  drawn  which  always  fall  in  a heavily  or 
lightly  shaded  region.  One  set  of  those  is  the  segmented  lines 
shown  in  Figure  11.  A surface  along  one  of  those  linos  is 
totally  dark  and  not  visible,  and  one  midway  between  these 
linos  is  relatively  bright.  However,  the  target  generally 
intersects  several  of  the  segmented  curves,  and  thus  appears 
alternately  dark  and  light. 

An  observer  viewing  the  master  grid  from  the  conver- 
gence point  sees  the  combination  of  the  master  grid  and  the 
superimposed  image  of  the  target  grid.  The  result  is  a moire 
fringe  pattern  which  changes  as  the  target  undergoes  out-of- 
plane motion.  In-plane  motion  of  the  target  does  not  affect 
the  target  grid  and  therefore  does  not  produce  a change  in  the 
fringe  pattern. 

The  straight  line  segments  drawn  on  the  moire  model  in 
Figure  11  connect  junction  points  (opposite  corners  of  heavily 
and  lightly  shaded  regions)  and  extend  ultimately  through  both 
the  source  point  and  the  convergence  point.  It  can  be  shown 
that  the  junction  points  for  one  segmented  line  lie  on  a smooth 
curve.  Data  reduction  was  greatly  simplified  by  the  deduction 
that  these  curved  lines  are  sections  of  ellipses.  This  deduc- 
tion is  explained  later. 

Figure  11  shows  the  scheme  for  relating  points  on  the 
curves  to  points  on  the  grids.  Arbitrary  reference  points  are 
chosen  at  and  X2  on  the  master  and  specimen  grids.  Rays 
and  A^  pass  through  the  grids  at  these  points  and  intersect 
on  the  i ellipse.  Rays  and  B2  pass  through  the  grids  at 
and  ^2’  is  observed  that  when  equals  X2Y2,  the 

rays  B^  and  B2  intersect  on  the  i^^  ellipse. 

Ray  B'2  passes  through  the  specimen  grid  at  a distance 
n/p  from  Y_  (n=l  is  illustrated).  Ray  B'^  intersects  ray  B. 
on  the  (i+n)^  ellipse. 
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The  curves  illustrated  in  Figure  11  are  called  dark 
curves  since  their  intersections  with  surfaces  produce  dark 
fringes.  As  the  dark  fringes  are  the  most  convenient  for 
measurement  and  the  dark  curves  are  the  simplest  to  illustrate « 
they  are  discussed  here  exclusively.  However,  the  analysis  can 
be  generalized  to  any  ellipse,  dark,  light,  or  some  fixed 
intermediate  shade. 

The  analycic  description  of  the  dark  curves  is  based 
on  the  moire  triangle  with  the  equivalent  specimen  and  master 
grids  drawn  normal  to  the  optical  axes.  To  show  that  the 
junction  points,  referred  to  above,  lie  along  ellipses  the 
moire  construction  shown  in  Figure  12  is  used.  The  moire 
triangle  and  the  equivalent  grids  are  shown  as  bold  lines. 


Figure  12.  Moire  Analytical  Diagram. 


The  triangle  is  described  by  three  parameters;  the  angle  between 
the  observation  axis  and  the  baseline,  a , and  the  angle  between 
the  projection  axis  and  the  baseline,  a , and  the  length  of  the 
baseline,  2L. 

The  location  of  the  dark  curves  and  their  characteris- 
tics may  be  formulated  as  follows.  An  arbitrary  ray,  A2,  is 
chosen  (shown  as  a dashed  line  in  Figure  12)  to  pass  through 
the  equivalent  specimen  grid  at  a distance,  g,  from  the  pro- 
jection axis.  This  ray  intersects  a complementary  ray, 
which  passes  through  the  equivalent  master  grid  also  at  a dis- 
tance g from  the  observation  axis.  The  distance  g can  be 
chosen  such  that  the  two  rays  intersect  on  one  of  the  dark 
curves  shown  in  Figure  11.  The  angles  between  these  rays  and 
the  baseline  are  B2  ^^d  Uj  as  shown  in  Figure  12.  A second 
pair  of  rays  B2  and  B^,  which  intersect  on  this  same  curve,  is 
illustrated  by  the  dotted  lines.  For  these  rays  to  be  on  the 
same  curve,  they  must  cross  both  grids  at  the  same  distance,  d, 
from  the  first  rays  as  demonstrated  previously  from  Figure  11. 
With  respect  to  an  origin  at  the  center  of  the  baseline,  the 
intersection  of  B2  and  B^^  is  the  point  (x,y)  . From  the  geo- 
metry, the  distance  y can  be  written  as 
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2L 


tan  y • tan  v 

I i 

tan  y + tan  y 
1 


(4) 


It  is  necessary  to  make  the  proper  substitutions  in 
Equation  4 so  that  y can  be  written  as  a function  of  x and  the 
physically  measurable  moire  parameters;  L,  a,  g,  f ^ , and 

f2*  The  first  step  is  to  make  the  substitutions 


'*1  = ‘M  “ 

y o ■=  « 2 + ^^2 


(5) 


where  aj  and  1x2  are  the  base  angles  of  the  moire  triangle  and 
the  angles  0i  and  Q2  are  taken  as  positive  if  they  represent 
clockwise  rotations  from  the  optical  axis.  Next,  Oy  and  O2 


must  be  replaced.  This  can  be  done  by  writing  O2  as  a function 
of  Oj  and  then  solving  0i.  Thus,  02  can  be  written  in  terms  of 
01  by  solving  the  equations 

tan  ©1  “ 

^1 

(6) 

tan  O2  * 

^2 


simultaneously  to  get 

tan  02 


2g  + tan  ©I 


f 


2 
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And  01  can  be  determined  from 


tan 


(ai  - Oi) 


-iL- 

L+x 


(8) 


After  expanding  tan  (ai  - Oi)  as  (tan  ai  - tan  Oi)/ 
(1  + tan  vx  1 tan  ©i),  Equation  (8)  can  be  written  as 


tan  0 1 


(L4-X)  tan  ui  - y 
y tan  ui  + (l+xT 


Finally,  Equation  (4)  can  be  rewritten  and  simplified  as 

Dy"+£y  (l+x)  +P\l+x)‘ 
where  the  coefficients  are  defined  as: 


(9) 


(10) 


A - ^ 9 tan  ai+  f^  tan  ix;^  tan  aj-  f^ 

B tan  «:'+  2g 

C = 0 (H) 

D = f,  tan  u;;+  f2  tan  aj  -2g  tan  vi  1 tan  a; 

£ = 2g  (tan  aj-  tan  01  o)  + ^ 

F = B 
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Equation  (10)  can  be  rewritten  in  the  form  of  the  general  ellipse; 


+ iSxy  + Qy^  + £y  + ^ = 

where  the  coefficients  are  defined  as 

^ tan  + £2  ^2 

jS  ^ 2q  (tan  - tan  a2)  + 

£ - ^2.  ^ ^2  *^1  ” “1  '^2 

a-  0 

£ = L (1-tan  tan  CX2)  -29  L(tan  + 

^ . L^/l 


0 (12) 


tan  a2) 


(13) 

tan  a.,  ) 


Equation  (12)  shows  that  the  curves  in  Figure  11  are 
ellipses,  and  the  coefficients  in  Equation  (13)  relate  the 
ellipses  to  the  moir€  parameters. 

Finge  data  are  reduced  by  measuring  the  parameters 
^1'  ^2'  ^*1'  ^^2  ' '^^lues  of  g,  separated  by  1/p,  are  used 

with  Equation (8 ) to  calculate  the  six  coefficients  of  any 
desired  ellipse.  The  displacement  of  a point  on  a moving  panel 
through  which  a given  fringe  passes  can  be  easily  calculated 
using  Equation  (12)  if  the  value  of  g is  known  and  the  x-position 
of  the  fringe  is  measured.  Equation  (12)  is  then  solved  for  y, 
the  distance  of  the  fringe  from  the  baseline. 

An  example  of  the  ellipse  generated  by  the  raoir^  data 
is  shown  in  Figure  13.  These  ellipses  were  based  on  the 
following  data:  f^^  » 6.51  cm,  f2  = 7.34  cm,  = 72'^,  a2  = 62.9°, 
L = 66.4  cm,  and  p = 78.7  lines/cm.  Fifty  values  of  g were 
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Figure  13.  Moire  Ellipses. 

used  ranging  from  0.127  to  0.495  in  increments  of  1/p.  The 
ellipses  were  numbered  from  1 to  SO.  Three  of  the  ellipses 
are  illustrated;  numbers  I#  25 « and  50.  Figure  14  illustrates 
the  tank  panel  shapes  calculated  for  shot  FTSB  using  this  method. 

3.2.1  Code  Description 

A computer  code  was  written  which  uses  the  theory  of 
moire  ellipses  to  calculate  out-of-plane  displacements  of 
panels.  As  described  above,  the  theory  of  moire  ellipses 
relates  the  fringe  positions  on  the  panel  to  intersections  of 
the  panel  with  a set  of  elliptic  cylinders.  These  ellipses  are 
completely  determined  by  six  parameters;  p,  the  pitch  of  the 
Ronchi  rulings;  2L,  the  length  of  the  base  of  the  moire  triangle; 
Oil  and  a2>  the  base  angles  of  the  moire  triangle;  and  f^  and 
f2»  the  distances  along  the  optical  axes  from  the  base  vertices  to 
the  equivalent  grids.  The  inputs  required  to  generate  the  ellipse 
information  and  the  outputs  from  the  code  are  described  below. 

Two  types  of  input  are  required  for  this  program. 
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Figure  14-  Panel  Shapes  Calculated  ior  Shot  FT5B  at  Various 
Times  After  Impact.  Interframe  time  is  0.165  ms; 
frame  1 occurs  between  t=0  and  t=0.165  ms. 

The  first  type  consists  of  measured  moir6  triangle  parameters; 
the  pitch  of  the  Ronchi  rulings,  the  base  length,  and  the  base 
angles.  The  remaining  required  information  are  the  lengths  f^ 
and  f p . These  distances  are  difficult  to  measure  but  can  be 
determin<jd  precisely  from  fiducial  pictures.  The  fiducial  pic- 
tures are  used  to  determine  the  coefficients  of  tv.'O  ellipses 
(called  calibration  ellipses) . Thus,  the  second  type  of  input 
consists  of  the  coefficients  of  the  calibration  ellipses,  and, 
from  these,  f,  and  can  he  determined.  Once  the  moir6  parameters 
are  known,  the  coefficients  of  any  desired  ellipse  can  be 
calculated . 

The  fiducial  displacement  data  are  determined  by  fixing 
a fiducial  to  the  target  v/ail  and  deterroini.ng  the  intersections 
of  the  moire  ellipses  with  the  fiducial.  An  illustration  of 
the  fiducial  is  shown  in  Figure  15.  Data  are  required  for  two 
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Figure  15.  Cross  Section  of  Fiducial  Step  Structure. 

adjacent  ellipses  which  intersect  the  fiducial  (the  calibration 
ellipses) . Each  set  of  calibration  ellipse  data  is  fit  to  the 
equation  for  an  ellipse.  Adjacent  ellipses  are  chosen  for  the 
calibration  data  and  are  arbitrarily  numbered  as  10  and  11. 

The  nvunbering  system  for  the  ellipses  must  meet  the  criteria  : 

1)  the  depth  covered  by  the  ellipses  must  be  enough  for  any 
anticipated  wall  displacement  (50  ellipses  are  currently  used) 
and,  2)  no  ellipse  can  be  labeled  0 or  less. 

An  iterative  scheme  is  used  to  determine  the  para- 
meters f^  and  The  distance  between  ellipses  is  determined 

from  the  calibration  ellipse  equations.  The  code  then 
determines  the  f^  and  f^  which  give  the  n ^asured  separation 
between  the  ellipses. 

The  ellipse  progreun  consists  of  the  main  program,  10 
subroutines,  and  three  function  subprograms.  The  main  program 
performs  no  calculations  but  calls  the  appropriate  subroutines. 
The  subroutines  are  described  briefly  below. 


1.  READIN 

All  data  is  input  through  the  same  subroutine. 

There  are  two  entry  points  to  this  subroutine  READANG  and 
READCAL . 

2 . CROSS  1 

This  routine  solves  the  equation  of  the  twc 
calibration  ollipsi>s  simultaneously  to  find  their  intersections. 
The  result  is,  in  general,  a fourth  degree  equation  in  y whose 
zeros  correspond  to  the  y coordinates  of  the  intersections  of 
the  ellipses.  CROSS  1 calls  ZEROS  to  calculate  these  inter- 
section points,  then  chooses  the  desired  pair. 

3 . ZEROS 

This  subroutine  solves  a fourth  degree  equation  in 
y giving  the  intersection  points  of  the  two  ellipses.  In 
general,  there  are  four  intersection  points  and  the  two  desired 
points  are  selected  in  CROSS  1. 

4 . CROSS  2 

In  the  event  that  the  two  ellipses  have  no  rotation 
or  translation  with  respect  to  the  baseline  (their  coefficients 
13  and  D are  both  zero) , the  simultaneous  solution  of  the  two 
calibration  ellipse  equations  yields  a quadratic  equation. 

CROSS  2 solves  this  quadratic  equation. 

5.  NEWORGN 

The  line  which  connects  the  intersection  of  the 
two  ellipses  is  defined  as  the  baseline  for  the  moire  triangle. 
The  midpoint  of  the  baseline  is  taken  as  the  origin  of  the 
moire  system.  The  coefficients  of  the  calibration  ellipses 
must  be  translated  and  rotated  to  align  with  the  baseline 
at  the  origin.  It  is  only  from  this  reference  system  that  the 
moire  parameters  can  be  used  to  generate  the  desired  ellipse.s. 
NEWORGN  locates  this  origin  and  executes  the  translation  and 
rotations. 

6 . PARAM 

This  subroutine  determines  the  moird  parameters 
needed  to  generate  the  desired  elipses.  One  of  its  principal 
functions  is  to  control  the  iteration  steps  in  the  search  for 


29 


values  of  and  f2  which  give  the  desired  separation  between 
ellipses. 

7 . GENCOF 

This  subroutine  calculates  the  coefficients  for 
the  desired  ellipse  based  on  the  moire  parameters.  The  ellipses 
are  numbered  from  1 to  50  with  the  largest  ellipse  being  number 

1. 

8 . REORGN 

The  coefficients  for  the  50  ellipses  are  trans- 
lated and  rotated  back  to  the  original  reference  system  of  the 
tank  wall  by  this  subroutine. 

9.  DISPLAC 

The  displacement  of  the  tank  wall,  y,  at  a fringe 
position  is  calculated  by  DISPLAC  which  solves  the  equations  of 
the  ellipses,  given  the  number  of  the  ellipse  and  x'  position 
of  the  fringe. 

10.  DYCALC 

This  subroutine  calculates  the  distance  between 
the  calibration  ellipses.  This  is  used  during  the  iteration 
procedure  used  to  find  the  correct  values  of  fj^  and  f2. 

3.3  EMPIRICAL  METHODS  OF  MOIRE  FRINGE  DATA  REDUCTION 

Two  methods  were  developed  for  empirically  converting 
fringe  shift  data  to  panel  displacement  information.  Both 
are  based  upon  study  of  the  fringes  produced  on  a fiducial 
surface.  The  first  uses  a step  surface  and  the  second  a wedge 
surface. 

3.3.1  Step  Method 

The  step  method  of  reducing  fringe  data  was  found  to 
be  the  most  useful  empirical  approach.  In  order  to  apply  this 
technique  a photograph  is  made  of  the  fringes  produced  on  a 
reference  surface  consisting  of  many  steps.  Figure  15  is  a 
diagram  of  the  fiducial  step  structure,  and  Figure  16  illus- 
trates the  fringes  on  this  fiducial  structure  for  shots  FT5B 
and  FT26A. 


Figure  16.  Fringes  Produced  on  Fiducial  Step  Structure  for 
Shots  FT5B  and  FT26A. 

Each  fringe  in  Figure  16  is  given  a designation,  and 
the  positions  on  the  photographic  record  at  which  each  fringe 
appears  on  each  step  are  measured.  The  positions  of  the 
troughs  of  the  fringes  are  determined  from  the  ramp  portion  of 
the  fiducial  structure.  In  this  way  a fringe"  map"  is  prepared. 
Figure  17  is  a demagnified  reproduction  of  the  fringe  map  pre- 
pared for  Shots  FT5B  and  FT26A.  In  the  figure,  values  of 
fringe  position,  x',  measured  on  the  film  along  a line  which 
was  a true  horizontal,  are  plotted  against  height,  h,  for  each 
fringe.  In  Figure  17  the  data  obtained  from  the  fiducial 
structure  have  been  extrapolated  to  obtain  fringe  location 
beyond  the  fiducial  structure.  The  extrapolation  is  facili- 
tated by  the  fact  that  within  the  accuracy  to  which  the  film 
can  be  read,  the  difference  in  height  between  any  two  fringes, 

A h,  is  constant  at  a given  value  of  x'.  In  other  words,  the 
fringes  are  all  congruent.  The  error  in  any  value  of  h in 
Figure  17  was  about  one-tenth  of  a fringe  separation,*  thus,  the 
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Figure  17.  Fringe  Map  Prepared  from  Data  of  Figure  16. 

accuracy  was  taken  to  be  ±0.2  mm  (0*008  in  } . 

Data  tables  were  prepared  for  x'  as  a function  of 
frame  number  for  each  fringe.  By  use  of  the  fringe  map,  these 
data  were  converted  to  h values  for  each  fringe  in  each  frame. 

Thus,  the  height  of  the  wall  surface  at  fringe  locations  in 
each  frame  was  determined,  with  the  aid  of  photographs  of 
a reference  grid,  x'  coordinates  were  transformed  to  x,  the 
horizontal  distance  parallel  to  the  undeformed  target  surface. 

These  data  are  plotted  and  fit  with  smooth  curves.  Such  a plot 
for  Shot  FT5B  data  is  shown  in  Figure  18.  The  fitted  curves 
are  drawn  dashed  in  the  region  of  the  center  of  the  fringe 
pattern,  where  the  fringes  are  farther  apart  and  the  panel 
surface  less  well  determined.  Because  the  data  points  are 
fitted  with  a smooth  curve,  it  is  possible  that  wall  deflections 
of  wavelength  shorter  than  ~20  mm  will  be  filtered  out  in  this  step 

The  error  in  initial  wall  shape  in  Figure  18  is  likely 
to  be  significant,  as  the  undeflected  target  panel  may  have 


been  at  an  angle  to  the  reference  surfaces  of  Figure  15.  How- 
ever a few  degrees  misalignment  will  have  an  insignificant 
effect  on  relative  displacement  measurements.  The  data 
reported  in  Section  V were  obtained  by  subtracting  the 
pre-impact  curve  from  the  various  post-ir»  "^ct  curves  in  Figure 
18. 

3.3.2  Wedge  Method 

Another  method  of  data  reduction  is  based  on  the  ob- 
servation that  the  slope  of  any  fringe  on  a non-vertical 
fiducial  plane,  which  intersects  the  undeformed  target  plane 
along  a horizontal  line,  is  a linear  function  of  x.  Thus,  we 
can  write 


X = +b 

where  C is  the  slope  of  the  fringe  (referenced  to  the  vertical) 
and  k and  b are  determined  from  a fiducial  photograph.  11  the 
angle  between  the  fiducial  plane  and  the  undeformed  target 


Figure  18.  Panel  Shapes  Derived  for  Fringe  Data  from  Shot 
FT5B  and  Fringe  Map  of  Figure  17.  Interframe 
time  was  165us. 


33 


plane  is  <{),  it  is  straightforward  to  show  that  for  any  one 
fringe 


dh  = tan  (j)  tan  ^ 
dx 


where  h and  x lie  along  that  fringe.  Integrating  gives 


h. . ,-h. 
1+1  1 


= m tan  4(  iln 


oos<=‘i-‘'> 


cos^^i+l"*^^ 


(14) 


the  displacement  increment  for  a fringe  which  moves  from 
x^  to  between  ith  and  i+lth  frames.  Since  Equation  (14)  only 

gives  a fringe's  displacement  relative  to  its  initial  position,  no 
useful  information  can  be  obtained  until  the  initial  tank  sur- 
face is  determined. 

This  technique  has  the  advantage  over  the  step  method 
that  data  reduction  involves  more  arithmetic  and  less  graphi- 
cal work.  Unfortunately;  however,  the  procedure  for  deducing 
the  initial  panel  shape  is  rather  imprecise,  for  it  requires 
adding  incremental  displacements  away  from  some  reference 
point  (such  as  a bolt) . If  the  panel  surface  is  well  surveyed 
just  prior  to  the  shot,  and  that  data  is  folded  into  the 
analysis,  the  wedge  method  would  probably  be  as  useful  as  the 
step  method. 
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SECTION  IV 

EXPERIMENTS  AND  RESULTS 


4 . 1 INTRODUCTION 

A series  of  impacts  into  water -filled  tanks  was  carried 
out  to  obtain  data  relevant  to  the  objectives  of  the  program. 
Impact  velocities  were  varied  from  1.18  to  2.38  km/s.  Out-of- 
plane displacement  data  were  obtained  for  bare  and  foam-backed 
entrance  and  side  panels.  Projectile  trajectory  and  shock  pro- 
pagation data  were  obtained  and  used  as  input  for  the  computa- 
tional treatments.  The  shots  for  which  data  was  obtained  are 
listed  in  Table  1.  Three  technical  areas  received  special 
attention  during  the  shot  sequence.  These  are  discussed  below. 

4.1.1  Projectile  Annealing 

The  projectiles  used  were  commercial  ball  bearings,  as 
described  in  Section  2.1.1.  The  initial  shots  in  the  program, 
FT1-FT3,  were  carried  out  with  the  bearings  in  the  as-received 
condition.  However,  it  was  found  that  these  bearings  broke  up 
on  impact  and  fragmentation  of  the  projectile  in  the  tank 
greatly  complicated  analysis.  In  order  to  eliminate  this  pro- 
blem, the  bearings  were  annealed.  They  were  annealed  by  heating 
to  400*C  and  holding  for  three  hours  followed  by  a one  hour 
cool  to  room  temperature.  A vacuum  furnace  was  used  for  the 
annealing  process. 

4.1.2  High  Velocity  Techniques 

Shots  FTlO-25  were  conducted  for  the  purpose  of 
developing  launch  techniques  for  high  velocity  projectiles.  It 
was  desired  to  achieve  velocities  in  excess  of  2.5  km/s.  How- 
ever, the  main  range  was  able  to  accommodate  barrels  up  to  only 
2 m in  length.  Within  the  constraints  on  barrel  length  imposed 
by  the  range,  such  high  velocities  were  found  to  be  impossible 
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except  vith  one-piece  lexan  sabots.  The  one-piece  sabots  could 
not  be  sufficiently  retarded  in  the  blast  tank,  and  impacted  the 
target  immediately  after  the  ball  bearing.  Double  impacts  such 
as  these  were  deemed  not  useful  for  the  progrenn.  However, 
techniques  were  developed  for  launch  of  14.3  mm  spheres  to 
velocities  in  excess  of  2.2  km/s.  The  kinetic  energy  of  such  a 
projectile  is  29  kj.  These  high  energy  impacts  were  employed 
for  study  of  side  panel  motion  (Shots  FT29-34) . 

4.1.3  Range  Instrumentation 

The  range  instrumentation  performed  generally  satis- 
factorily during  the  program.  The  flash  x-ray  records  were  found 
to  be  extremely  valuable  for  diagnosing  projectile  performance 
during  sabot-development  teat  shots.  A typi^cal  radiograph 
of  the  projectile  in  flight  is  shown  in  Figure  19.  This  shows 
a sphere  as  it  just  touches  a trigger  screen.  The  sphere  is  not 
broken  and  all  sabot  parts  have  been  satisfactorily  stripped 
away. 


4.2  DESCRIPTION  OF  PANEL  DAMAGE 

4.2.1  7075-T6  Aluminum  Entrance  Panels 

In  Shots  FT2-4,  the  entrance  panels  were  made  from 
7075-T6  aluminum  alloy.  The  damage  done  to  these  panels  by 
hydrodynamic  ram  is  shown  in  Figures  20-22 . Cracks  initiated 
at  the  entrance  site  rapidly  and  propagated  to  the  edges  of  the 
plates.  In  Shot  FT3  one  crack  traversed  the  field  of  view  of 
the  framing  camera  recording  the  moire  fringe  pattern.  The 
average  crack  velocity  over  the  first  450  mm  of  motion  (to  the 
fasteners)  was  310  ± 40  m/s.  This  is  about  1/5  of  the  sound 
velocity  for  water. 

4.2.2  2024-T3  Aluminum  Entrance  Panels 

A major  purpose  of  the  program  was  to  obtain  displace- 
ment data  for  panels.  Thus,  it  was  necessary  to  measure  panel 
motion  at  several  times  before  failure  or  measure  motion  of 
panels  that  did  not  fail.  It  was  also  desired  to  remain  in  the 
impact  range  characteristic  of  high-velocity  fragments. 
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Figure  19. 


Plash  X-Radiograph  of  11. i nun  Diameter  Projectile 
sLt  PT26a!''^®  Second  Mylar  Foil  Switch  in 


Figure  20.  Panel  from  Shot  PT2  after  Impact.  Panel  ia  .1.6  mm 
thick  7075-T6  aluminum.  Impact  was  11.1  mm 
diameter  sphere  at  1.73  km/s. 
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Figure  21.  Panel  from  shot  FT3  after  Impact.  Panel  is  1.6  nmi 
thick  7075-T6  aluminum.  Impact  was  11.1  n\m  dia- 
meter sphere  at  1.73  knv's.  Scale  divisions  on  panel 
are  51  n\m  apart. 


Figure  22.  Panel  from  Shot  FT4  after  Impact.  Panel  is  1.6  nun 
thick  7075-T6  aluminum.  Impact  was  14.3  nun  dia- 
meter sphere  at  1.61  knv's.  Scale  divisions  on 
panel  are  51  ;'’m  apart. 
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(e.g.,  1.5  kiti/s)  . Failures  were  observed  to  occur  very  early 
with  7075-T6  entrance  panels.  Since  2024-T3  aluminum  alloy  is 
substantially  more  ductile  than  7075-T6  alloy,  entrance  panels 
made  from  this  material  were  likely  to  be  more  resistant  to 
massive  cracking.  This  turned  out  to  be  the  case.  Therefore, 
2024-T3  entrance  panels  were  used  for  the  remainder  of  the  program. 

It  was  found  that  11.1  mm  diameter  spheres  at  velocities 
less  than  about  1.6  km/s  did  not  produce  any  cracks  in  the  en- 
trance panels.  Cracks  produced  by  higher  velocity  shots  against 
2024-T3  panels  were  relatively  short.  Figures  23  and  24  show  ex- 
amples of  damage  produced  in  2024- T3  alloy  entrance  panels.  A 
sequence  of  photographs  which  illustrate  damage  caused  by  in- 
creasing projectile  velocities  with  the  28  £ tank  is  shown  in 
Figure  25  through  30.  At  1.21  and  1.60  km/s  impact  velocity  no 
damage  resulted  (other  than  the  round  entrance  perforation) . At 
1.94  km/s  small  cracks  formed  at  the  impact  site,  but  did  not  grow. 
At  2.21  km/s  four  cracks  grew  which  nearly  cut  the  panel  into  four 
pieces.  At  2.23  km/s  the  panel  was  cut  into  four  pieces.  However, 


Figure  23.  Panel  from  Shot  FT5  after  Impact.  Panel  is  1.6  mm 

thick  2024-T3  aluminum  Impact  was  11.1  mm  diameter 
sphere  at  1,53  km/s.  Scale  division  on  panel  are 
51  mm  apart. 
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Figure  28.  Panel  from  Shot  FTA6  , 1.6  mm  thick  2024-T3  aluminum. 

Impact  was  11.1  mm  diameter  sphere  at  2.21  km/s . 
Circular  bend  was  Caused  by  Collision  with  Kantje 
Hardware  after  Impact. 


Figure  29.  Panel  from  Sliot  FTAl  0 , 1.6  nun  thick  2024-'r3  aluminum. 

Impact  was  11.1  mm  diameter  sphere  at  2.23  km's. 
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Figure  30.  Panel  from  Shot  FTA12,  1.6  nvm  thick  2024-T3  aluminum 
backed  with  11  ram  AVCO  ballistic  foam.  Impact  was 
11.1  mm  diameter  sphere  at  2.38  km/s. 

at  2.38  km/s  with  11  mm  of  ballistic  foam,  no  cracking  occurred. 
This  illustrates  dramatically  the  ability  of  foam  backing  to  defeat 
hydrodynamic  ram. 

In  summary,  after  a hydrodynamic  ram  event,  2024-T3  alloy 
entrance  panels  generally  are  bulged  outward  around  the  entrance 
hole.  Figure  31  shows  typical  profiles.  The  tensile  strains  asso- 
ciated with  the  enlargement  of  the  entrance  hole  radius  may  give 
rise  to  cracks  which  then  propagate  due  to  the  loading  applied 
during  the  drag  phase. 

4.2.3  Effect  of  Foam  Protection  on  Hydrodynamic  Ram  Damage 

Previous  investigations  showed  that  hydrodynamic  ram  damage 

to  fuel  cells  can  be  significantly  reduced  by  the  presence  of  en- 
. (17-19) 

ergy-absorbing  substances'  . The  most  promising  technique 

appears  to  be  a semi-rigid  foam  interface  between  the  wall  and  the 
fluid.  The  usefulness  of  two  types  of  foams  for  defeat  of  hydro- 
dynamic  ram  was  investigated  in  the  present  program:  AVCO  AX  5052- 
2.5  Thermarest  ballistic  foam  and  styrofoam. 


SHOT  FT  5 

11.1  mm  DIA. 
1.53  km/s 


SHOT  FT5B 

11.1  mm  DIA. 
1.55  km/s 


SHOT  FT  26A 

11.1  mm  DIA. 

I. 52  km/s 
38  mm  FOAM 

SHOT  FT  28 

II. 1  mm  DIA. 
2.18  km  /s 
33  mm  FOAM 

SHOT  FT  8 
14.3  mm  DIA. 
1.59  km/s 
38  mm  FOAM 


SHOT  FTA4 

11.1  mm  DIA 
1.60  km/s 


Figure  31.  Profiles  Through  Impact  Hole  of  Selected 
2024-T3  Alloy  Entrance  Panels. 
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Most  tests  of  foam  were  performed  with  the  AVCO  balli- 
stic foam.  In  Shots  FT7,  FT8,  FT26,  FT26A,  FTAll,  and  PTA12, 
varying  thicknesses  of  AVCO  foam  were  bonded  to  the  entrance 
panel  with  epoxy.  Initially  76  mm  thick  foam  was  employed, 
and  that  thickness  proved  extremely  effective.  The  foam 
thickness  was  reduced  in  successive  shots  to  determine  a mini- 
mum useful  foam  thickness.  In  the  final  shots,  the  foam 
thickness  was  only  11  mm  and  yet  the  hydrodynamic  ram  mechanism 
was  completely  defeated.  None  of  the  foam-backed  entrance  panels 
were  significantly  damaged  in  these  experiments.  Comparison  of 
Figure  30  with  Figure  29  provides  an  indication  of  the  effective- 
ness of  a thin  foam  layer  behind  the  entrance  panel . The  same 
projectile  which  destroyed  a panel  with  no  foam  backing  caused 
insignificant  damage  to  the  panel  after  it  was  protected  with  an 
11  mm  thickness  of  rigid  foam. 

The  effectiveness  of  styrofoam  was  evaluated  in  Shot 
FT9 . The  entrance  panel  was  7075-T6  alloy.  The  extensive 
cracking  which  occurred  on  bare  walls  was  completely  eliminated 
by  50  mm  of  styrofoam. 

Post  shot  observations  of  the  foam-backed  panels 
revealed  that  the  foam  was  typically  crushed  within  a radius 
of  120  to  150  mm  centered  on  the  impact  hole  in  the  panel. 

The  severely  damaged  region  was  only  10  to  30  iim\  in  radius  on 
the  side  of  the  foam  adjacent  to  the  fluid.  The  large  damage 
next  to  the  panel  probably  resulted  from  dissipation  of  the 
shock  wave  in  the  foam.  Fluid  shock  effects  are  reduced  when 
the  projectile  strikes  the  foam  backed  fluid  because  the  peak 
shock  pressure  is  lower  (see  Figure  1) , and  the  foam  permits 
rapid  release  of  the  shock  presures. 

Some  insight  into  the  mechanics  of  projectile,  foam, 
panel,  and  water  interactions  can  also  be  obtained  from 
analysis  of  Figure  31.  Note  that  the  permanent  deformation 
which  resulted  in  FT5,  FT5B,  and  FTA4  (which  represent  nearly 
identical  impact  configurations)  was  completely  absent  in 
FT26A  (the  same  impact  configuration  with  the  addition  of  foam) . 
The  relatively  long  wavelength  of  the  deformation  suggests 
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that  it  is  caused  by  drag  phase  pressure.  Therefore,  the  chief 
effect  of  the  foam  was  apparently  to  rigidly  shield  the  wall 
from  the  drag  phase  pressure  pulse.  As  the  ball  size  was 
increased  (FT8) , the  response  of  foam-backed  panels  did  not 
change.  In  the  high  velocity  Shot  FT28,  large  local  deformation 
occurred.  This  was  probably  due  to  the  shock  phase  pressure, 
which  increases  with  impact  velocity,  and  which  propagates  for 
a short  distance  into  the  foam.  However,  the  lack  of  a larger 
wavelength  bulge  suggests  that  the  overall  rigidity  of  the  foam 
shielded  the  entrance  panel  from  cavity  phase  pressure.  There 
are,  therefore,  two  mechanisms  by  which  the  rigid  foam  defeats 
hydrodynamic  ram  at  entrance  panels:  the,'  crushing  at  the  foam 
attenuates  the  shock  pressure;  and  the  flexural  rigidity  of  the 
foam  spreads  the  drag  phase  pressure. 

One  shot  was  carried  out  with  ballistic  foam  on  a side 
panel,  Shot  FT34 . No  side  panel  failures  or  large  plastic 
deformations  had  been  observed,  and  thus  it  was  felt  that  a 
closer  trajectory  would  be  required  to  evaluate  the  foam.  The 
trajectory  was  moved  horizontally  about  150  mm  closer  to  the 
observed  panel.  No  failure  occurred,  but  moire  fringe  data  were 
obtained.  The  displacement  results  are  reported  in  Section  4.3. 

Because  of  the  long  trajectory  in  the  water  and  low 
projectile  masses,  rear  panels  were  only  occasionally  impacted 
with  sufficient  force  to  cause  failure.  The  ability  of  foam 
to  reduce  rear  panel  damage  was  nevertheless  investigated  in 
Shots  FT31-FT34.  In  all  four  shots  the  rear  panel  was  cracked 
but  not  punctured.  Neither  50  ituri  of  styrofoam  nor  12  mm  of 
ballistic  foam  appeared  to  reduce  the  extent  of  cracking.  The 
difference  in  foam  effectiveness  for  entrance  and  side  panels 
may  be  due  to  the  fact  that  the  fluid  velocity  vector  is 
parallel  to  or  away  from  the  entrance  panel  but  perpendicular 
to  or  toward  the  exit  panel. 

4.3  MEASUREMENT  OF  WALL  DISPLACEMENT 

Out-of-plane  displacement  data  were  reduced  for  six 
shots,  as  of  the  writing  of  this  report.  The  shots  reduced 
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were  FT2,  FT5B,  FT6,  FT8,  FT26A  and  FT33.  FT33  was  a side  panel 
shot  but  the  remaining  five  were  entrance  panel  shots.  The  ''*iase'' 
shot  was  FT5B;  a 11.1  mm  dicuneter  sphere  impacted  at  1.53  km/s 
against  a 2024-T3  aluminum  panel.  FT2  was  similar  to  FT5/  except 
the  panel  was  7075-T6  aluminum.  FT26A  was  similar  to  FT5,  except 
the  panel  was  backed  by  38  mm  of  ballistic  foam.  FT6  was  the 
same  as  FT5,  except  that  a 14.3  mm  diameter  sphere  was  used. 

FT8  was  the  same  as  FT5  except  a 14.3  mm  diameter  sphere  was 
used  and  the  panel  was  backed  by  38  mm  of  ballistic  foam. 

4.3.1  Entrance  Panel  Measurements 

Figure  32  illustrates  8 successive  frames  from  the 
base  Shot  FT5B.  The  target  point  is  the  cross  on  the  tank 
surface  in  frame  "I.  Impact  occurs  between  the  two  frames 

labeled  -1  and  0.  The  interframe  time  is  164  ms.  In  frame  -1, 
displacement  of  the  center  of  the  panel  has  caused  the  center 
of  the  fringe  pattern  to  shift.  By  frame  2,  the  impact  flash  has 
dissipated  sufficiently  to  view  most  of  the  target.  Starting 
in  frame  3 and  continuing  through  frame  5 a discontinuity  in 
the  fringe  positions  appears  to  propagate  away  from  the  impact 
zone.  In  frames  5 to  7,  a second  discontinuity  appears  to  follov; 
the  first.  In  later  frames  reflections  of  these  motions  from 
the  bottom  of  the  tank  are  apparent. 

The  moire  fringe  records  for  this  shot  were  read  along 
horizontal  and  vertical  lines  through  the  impact  point.  The  data 
reduction  technique  used  was  that  described  in  Section  3.3.1. 
Figures  33  and  34  show  the  tank  displacement  determined  from 
the  moire  fringe  data.  As  discussed  in  Section  3.3,  the  tank 
displacement  given  in  the  figures  may  be  lacking  some  high 
spatial  frequences,  since  data  is  only  available  at  fringe 
positions  and,  during  the  data  reduction  procedures,  a smooth 
curve  was  dravm  through  the  measured  fringe  positions.  The 
accuracy  in  displacement,  6 (Ah),  was  generally  ± 0.4  mm.  The 
accuracy  was  somewhat  reduced  over  the  panel  region  on  which 
the  center  of  the  fringe  pattern  occurred  because  fringes  are 
farther  apart.  For  example,  for  100  s x s 200  mm,  6 (Ah)  is 
probably  about  0.5  mm.  The  effect  of  the  accuracy  in  Ah  is 
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Figure  33.  Out-Of-Plane  Displacement  (Measured  from  Pre-Impact 
Surface) . Error  bars  indicate  approximate  uncer- 
V,  tainty  in  displacement.  Time  uncertainty  is 

± 0.08  ms.  Data  are  for  entrance  panel  from  Shot 
FT5B,  an  11.1  mm  diameter  ball  striking  a 2024  T3 
panel  at  1.53  km/s. 
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that  the  curves  in  Figures  33  and  34  may  be  redrawn  within  bands 
of  ± 6 (Ah)  without  conflicting  with  the  moire  fringe  data. 

In  FT5B  the  maximum  displacement  at  the  impact  zone 
was  approximately  29  mm,  which  was  attained  at  about  1.7  ms  after 
impact.  Outward  movement  of  the  flanks  of  the  panel  continued 
beyond  this  time,  while  the  peak  displacement  began  to  fall. 
Deflections  on  the  line  below  the  impact  point  are  about  10% 
higher  than  at  comparable  distances  measured  horizontally  from 
the  impact  point.  Motion  near  the  panel  edges  starts  between 
200  and  400  us  after  impact.  The  "pucker"  of  the  panel  in  the 
entrance  hole  region  was  established  by  740  us  after  impact, 
but  was  absent  before  this  time.  The  initiation  of  major  motion 
of  the  target  panel  propagates  across  the  panel  with  a speed  of 
300  ± 30  m/s.  A smaller  negative  displacement  wave  appears  to 
propagate  to  about  100  mm  distance  at  about  1 ms  after  impact 
with  a speed  of  122  ±24  m/s.  (Displacement  front  motion  can  be 
measured  directly  from  the  moire  framing  camera  pictures) . 

Displacement  data  for  the  corresponding  foam  shot, 
FT26A,  are  presented  in  Figures  35  and  36. 

The  panel  motion  in  FT26A  was  smaller  in  this  shot. 

Peak  displacement  '“'as  only  about  22  mm  compared  to  29  mm  for 
FT5.  Initial  panel  movement  was  considerably  slower,  but  the 
maximum  displacement  in  the  impact  zone  was  still  reached  at 
about  1.7  ms  after  impact.  The  curvature  of  the  outward  moving 
panel  was  also  less  in  this  shot.  The  major  motion  discontinuity 
was  of  smaller  amplitude  than  that  of  FT5B,  and  only  propagated 
at  170  ±14  m/s. 

Figure  37  presents  the  displacement  data  for  a similar 
shot  with  a larger  diameter  sphere,  FT8.  The  profiles  are  very 
similar  to  those  of  Figure  36.  FT6  is  the  companion  large  pro- 
jectile shot  with  no  foam.  That  panel  failed,  as  is  shown  in 
Figure  24.  The  displacements  along  the  horizontal  scale  to  the 
left  of  the  impact  site  are  shown  in  Figure  38.  In  this  shot, 
impact  flash  from  a foil  switch  obscured  the  area  within  75  mm 
of  the  impact  site.  Beyond  that  region,  in  spite  of  the  failure 
which  occurred  elsewhere,  the  displacement  along  this  parti- 
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HORIZONTAL  DISTANCE  FROM  IMPACT  POINT  (mm) 


Figure  35.  Same  as  Figure  33,  But  for  Shot  FT26A,  an  11.1  mm 
Diameter  Sphere  Striking  a 2024-T3  Panel  Backed  by 
38  mm  of  Ballistic  Foam  at  1.52  km/s. 


VERTICAL  DISTANCE  FROM  IMPACT  POINT  (mm) 


Figure  36.  Same  as  Figure  35,  But  for  Vertical  Line  Through 
Impact  Site. 


Figure  37,  Same  as  Figure  33,  Except  for  Shot  FT8,  a 14.3  mm 
Diameter  Sphere  Striking  a 2024-T3  Panel  Backed  By 
38  mm  of  Ballistic  Foam  at  1.59  km/s. 


HORIZONTAL  DISTANCE  FROM  IMPACT  POINT  (mm) 


Figure  38.  Same  as  Figure  33,  Except  for  Shot  FT6,  a 14.3  mm 

Diameter  Sphere  Striking  a 2024-T3  Panel  at  1.52  km/ 


cular  line  is  qualitatively  similar  to  that  of  the  unfailed  panel 
of  FT5.  The  major  difference  between  these  two  shots  is  the 
relatively  smaller  amount  of  displacement  near  the  fasteners  of 
the  failed  panel. 

Figure  39  presents  entrance  panel  displacement  data 
from  FT2.  This  panel  was  made  of  7075-T6  aluminum;  a post- 
impact photograph  of  the  failed  plate  is  shown  in  Figure  20. 

For  early  times  near  the  entrance  region  the  movement  is  similar 
to  that  of  Figure  35.  However,  the  peak  displacement  at  late 
times  is  very  large,  due  to  petalling. 


Figure  39.  Same  as  Figure  33,  Except  for  Shot  FT2,  an  11.1  mm 

Diameter  Sphere  Striking  a 7075-T6  Panel  at  1.73  km/s. 
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4.3,2 


Side  Panel  Measurements 


Moir6  fringe  data  from  the  side-panel  Shot  FT33  were 
reduced  along  a hori:eontal  line  at  trajectory  height  by  the 
technique  described  in  Section  3.3.1,  and  along  a horizontal  lino 
150  mm  below  trajectory  by  the  technique  described  in  Section 
3.3.2.  The  out-of-plane  panel  displacement  at  various  times  are 
displayed  in  Figures  40  and  41.  The  origin  on  these  plots  is 
the  center  of  a fastening  bolt. 

The  uncertainty  in  wall  shape  in  Figures  40  and  41  is 
about  i 0.25  isur. . Thus,  the  inflection  points  in  the  central 
regions  of  the  displacement  curves  seem  to  bo  real.  The  dip  in 
the  curves  near  the  fasteners  is  believed  to  be  due  to  twisting 
of  the  edge  angle  bracket  oau^^ed  by  front  panel  loading. 
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Figure  40.  Sai  ■»  as  Figure  33,  but  for  Side  Panel  of  Shot  FT33, 
a 14.3  ram  Diameter  Sphere  Striking  at  1.83  km/s. 
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Figure  41.  Same  as  Figure  40,  Except  for  a Line  150  mm  Below 
Trajectory . 

The  peak  displacement  experienced  by  side  panels  is 
not  reached  until  relatively  late — about  7 ms  in  FT33  and  6 ms 
in  FT32.  This  is  long  after  the  projectile  has  exited  the  tank. 
The  maximum  displacement  observed  is  about  25  mm.  This  is 
comparable  to  the  entrance  panel  displacement  in  FT5B,  even 
though  the  energy  of  the  projectile  in  shot  FT33  was  almost 
three  times  that  of  the  projectile  in  shot  FT5B. 

4.4  PROJECTILE  DRAG  DETERMINATION 

Figures  42  and  43  show  sequences  of  frames  from  typical 
film  records  of  projectile  and  cavity  motion  in  fluid  filled 
tanks.  Data  vere  available  from  the  large  tank  shots,  (as 
shown  in  Figure  42}  and  also  from  the  shock  tank  shots  (as 
shown  in  Figure  43)  at  early  times.  These  data  were  used  to 
determine  the  projectile  velocity  and  the  cavity  expansion  and 
decay  rates.  The  projectile  velocity  was  used  to  determine 
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coefficient  of  drag,  a parameter  which  was  used  as  input  to 

the  fluid  drag  code. 

The  coefficient  of  drag  is  defined  by  the  relation 
du  _ pu^ 


-m 


dt 


2 ^ 


(15) 


where  m and  u are  the  projectile  mass  and  velocity,  p is  the 
density  of  the  fluid,  a is  the  projectile  cross-section  area  and 
Cj^  is  the  coefficient  of  drag.  The  drag  force  on  a body  is, 
therefore,  proportional  to  the  stagnation  pressure  and  the  cross- 
sectional  area  of  the  body  normal  to  the  projectile  velocity.  The 
coefficient  of  drag  is  a constant  and  in  general  depends  on  the 
Reynolds  number  and  mach  number,  as  shown  in  Figures  44  through 
46.  The  data  in  these  figures  were  obtained  for  projectiles  in 
air. 


4.4.1  Projectile  Reynolds  Numbers 

The  Reynolds  number  for  a sphere  is  given  by 

Re  = Pf""^ 


(16) 


where  p^  is  the  fluid  density,  D is  the  diameter  of  the  sphere, 


Figure  44.  Drag  Coefficient  vs.  Reynolds  number  for  Spheres 
(Taken  from  Prandtl^^). 
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Figure  45.  Drag  Coefficient  vs.  mach  number  for  Spheres  (Taken 
from  Shapiro^**). 


Figure  46.  Reynolds  Number  vs.  Velocity  for  Spheres. 
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The  Reynolds  numbers  for  11.1  and  14.3  mm  diameter  spheres  are 
plotted  in  Figure  46. 

Comparing  Figure  46  to  Figures  44  and  45  one  sees  that  the 
initial  drag  coefficient  of  the  projectiles  (which  are  slightly 
supersonic)  will  be  0.6  to  0.8.  By  the  time  a projectile  is  mid- 
way across  the  large  tank  (0.5  m travel)  its  velocity  has  usually 
dropped  to  a little  less  than  half  its  initial  value.  Thus,  Cj^ 
is  expected  to  fall  to  about  0.4.  The  drag  coefficient  will  then 
continue  to  decrease,  unless  the  velocity  falls  below  about  100  m/s. 
At  that  velocity,  the  critical  Reynolds  number  regime  is  reached 
and  Cjj  increases  sharply.  The  observed  values  of  are  con- 
sistent with  these  considerations. 


4.4.2  Drag  Measurements 

To  determine  the  coefficient  of  drag,  the  projectile 
position  and  time  data  were  least-squares  fit  to  the  equation 

^ &n  (icUit+1)  ,, 


where 


K = 


Pf  Cp 


(18) 


and  u^  is  the  initial  velocity  of  the  projectile  in  the  fluid. 

In  Equation  (18) , and  p are  the  fluid  and  projectile  densities, 
and  R is  the  projectile  radius.  Equation  (17)  follows  from  Equation 
(15)  after  two  successive  integrations. 

The  value  of  u^  was  determined  before  the  large  tank  data 
was  curve  fit  to  Equation  (17).  This  determination  was  based  on 
displacement  vs.  time  information,  obtained  from  the  shock  tank 
shots.  These  measurements  provided  numerous  data  points  for  the 
first  5 to  10  cm  of  projectile  travel  through  the  fluid.  The 
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projectile  velocity  was  found  to  be  essentially  constant  from 
about  10  to  50  ys  after  impact.  Plots  of  initial  velocity  vs.  the 
impact  velocity  are  shown  in  Figure  47.  For  small  spheres,  vel- 
ocity lost  during  impact  was  approximately  a linear  function  of 
the  impact  velocity  given  by 

u -u.  = 0.92  u -1.1  (19) 

or  o 

where  u is  expressed  in  km/s.  Alternatively,  Equation  (19)  can 
be  written 

u^  = 1.1  t 0.08  u^  (20) 

Only  one  large  sphere  data  point  was  available,  so  no  conclusion 
could  be  drawn  concerning  the  behavior  of  large  spheres. 

Data  were  least-squares  fit  to  Equation  (17)  using  a u^ 
calculated  from  Equation  (20)  . The  resulting  drag  coefficients 
were  used  in  the  calculations.  The  range  of  variation  was  small, 
however,  ranging  from  0.32  to  0.36  and  increasing  with  impact 
velocity.  Figure  48  Compares  an  actual  trajectory  with  the  pre- 
dicted one  for  a typical  case  (FT5) . 

Another  important  result  is  evident  in  Equation  (20) . 

For  small  spheres,  the  initial  velocity  of  the  projectile  is 
relatively  insensitive  to  variations  in  impact  velocity.  Thus, 
the  high  velocity  projectiles  surrender  increasing  amounts  of 
kinetic  energy  immediately  at  impact.  This  energy  is  absorbed 
in  the  strong  shocks  associated  with  high  velocity  projectile 
impacts . 

There  is  evidence  that  the  has  a value  of  approximately 
0.6  in  the  region  within  about  15  cm  of  the  entrance  panel.  Values 
in  this  range  were  computed  by  CRT  using  the  AFTON  code.  The 
shock  tank  data  also  imply  large  values  of  but  those  data  are 

considered  unreliable  for  acceleration  measurements  due  to  the 
"lensing"  effect  of  the  shock  wave. 

4.5  SHOCK  VELOCITY  MEASUREMENTS 

The  propagation  speeds  of  impact-induced  shock  waves  in 
the  fluid  were  measured  using  the  shook  tank.  Peak  pressure  can 
be  calculated  directly  from  these  data  by  using  the  Hugoniot 
relationship : 
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p = Pq  UU, 


(21) 


(where  u and  u are,  respectively,  the  particle  velocity  and  the 
shock  propagation  velocity,  measured  relative  to  the  material 
into  which  the  shock  is  moving)  and  the  known  Hugoniot  of  water 

U = 1.51  + 1.85  u (in  km/s) . (22) 

Values  of  peak  shock  pressure  calculated  from  (21)  and  (22)  and 
the  measured  values  of  U along  the  trajectory  for  Shot  FTA9  are 
shown  in  Figure  49. 


Figure  49.  Peak  Shock  Pressure  Calculated  from  Shock  Velocity 
Data  from  Shot  FTA9. 
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For  present  purposes,  the  shock  posit  ion- tijme  data  were 

(7) 

fit  to  an  equation  provided  by  Yurkovitch'  ; 

X * kt*^'  (23) 

where  x and  t are  shock  front  position  and  time,  and  k and  are 
parameters  derived  from  measurements.  The  properties  of 
Equation  (23)  and  the  procedure  used  for  deriving  k and  are 
discussed  in  Section  5. 

The  effect  of  pressure  relief  on  the  front  tank  wall  can 
be  seen  by  comparing  shock  propagation  (x,  t)  data  into  the 
fluid  with  that  along  the  wall.  For  example,  Figure  50  shows 
(x,  t)  data  from  Shot  FTA9.  The  shock  propagation  velocities 
in  the  two  propagation  directions,  at  30  mm  in  the  tank,  are 
1.87  ±0.04  and  1.54  ±0.07  km/s,  respectively.  The  corresponding 
pressures  are  37.6  and  26,3  MN/m^. 


Figure  50.  Shock  Position  vs.  Time  for  Shot  FT9. 
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The  degree  of  pressure  relief  near  the  entrance  panel 
is  larger  for  foam-backed  panels.  In  FTA12  (38  nun  foam-backed 
panel)  after  mm  travel  the  shock  speeds  into  the  fluid  and 
along  the  wall  are,  respectively,  1.78  ±0.09  and  1.27  ±0.25  km/s. 
Corresponding  pressures  are  21.3  and  ^-0  MN/m^. 

4.6  PRESSURE  MEASUREMENTS 

Pressure  transducers  were  included  in  most  of  the  shots 
in  order  to  provide  data  to  check  code-generated  pressure  values. 

The  pressure  data  were  recorded  on  wide  bandwidth  oscilloscopes. 

Time  of  arrival  was  usually  determined  from  a fiducial  mark 
superimposed  on  the  data  channel. 

In  many  cases  the  initial  pressure  spike  was  off-scale  on 
the  oscillograph.  In  addition,  high  frequency  noise  (-100  kHz)  from 
mechanical  resonances  was  superimposed  on  many  data  traces. 

Figure  51  illustrates  a typical  example,  taken  from  Shot  FT31. 

The  transducer  was  situated  near  the  center  of  the  tank.  The 
initial  arrival  is  off  scale,  but  the  trace  returns  about  30  ps 
later.  There  is  a pressure  plateau  of  about  1.7  MN/m*,  and  the 
total  duration  of  the  pulse  is  about  150  ps.  This  transducer 
underwent  a slight  zero  shift,  which  occurred  in  about  one-third 
of  the  records.  At  later  times  there  is  a 2 kHz  oscillation. 

This  also  occurred  in  many  records,  and  is  a little  too  fast  to 
be  ascribed  to  standing  pressure  waves  in  the  fluid.  More 
likely,  it  is  caused  by  a resonance  in  the  transducer-support 
apparatus. 

For  many  records,  the  high  frequency  components  were  removed 

from  the  transducer  output  with  30  kHz  cut-off  filters.  Two  typical 

filtered  records  are  shown  in  Figure  52.  The  initial  peaks  may 

have  been  slightly  clipped,  since  their  rise  times  are  '^>10  ps.  The 

remainder  of  the  traces  are  a great  deal  easier  to  read  than  that 

of  Figure  51.  Prom  the  lower  trace  of  Figure  52,  which  is  from  the 

transducer  nearest  the  rear  panel  in  Shot  FT34/  it  can  be  seen  that 

2 

the  initial  pressure  spike  has  an  amplitude  of  v6.5  MN/m  and  a 
duration  of  '\>100  ps.  The  pressure  then  falls  to  zero  over  300  ps.  Of 
the  total  impulse  contained  in  the  pulse,  about  half  is  in  the  leading 
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Figure  51.  Shot  FT31,  Transducer  T3  (431  nun  from  Entrance  Panel, 
100  nun  Above  Trajectory)  . Scope  is  triggered'  at 
Impact.  0.49  ms/div  sweep  rate.  1.4  MN/m"/div  ver- 
tical scale. 


Figure  52.  Shot  FT34  Upper  beam:  Transducer  T4  (7  37  nun  from^ 

Entrance  Panel,  180  nm)  From  Trajectory)  1.35  MN/m'Vdiv 
Lower  Beam:  Transducer  T5  (983  nun  from  Entrance  Panel, 
180  mm  From  Trajectory),  1.6  MN/m^/div,  1.5  ms  de- 
layed sweep.  0.49  ms/div  sweep  rate. 
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spike.  This  indicates  that  even  deep  within  the  tank,  hydrody- 
namic ram  loading  is  a very  rapid  event  measured  in  tenths  of 
milliseconds. 

Additional  discussion  of  pressure  records  for  individual 
experiments  ca  be  found  in  Section  V.  There  the  measured 
pressures  are  compared  with  those  predicted  by  numerical  techniques. 


SECTION  V 
CALCULATIONS 


Two  approaches  have  been  taken  to  computer  prediction  of 
fuel  tank  wall  motion.  The  first  approach  was  based  on  a finite 
element  computer  code,  BR-lA,  originally  developed  by  Northrop 
Corpoxation  to  calculate  tank  wall  response  to  an  internal 
blast.  This  code  was  modified  to  account  for  fluid  loading  of 
the  tank  walls,  BR-lA  (HR).  A principal  input  is  the  fluid 
pressure  that  would  exist  at  < ho  boundary  if  the  walls  were 
absent.  These  pressures  were  generated  from  two  computer 

routines;  the  description  of  the  pressures  in  the  fluid  due  to 

, (14) 

projectile  drag  uses  a code  developed  by  Lundstrom  (fluid 

drag  code)  and  the  uressure  profile  behind  iJu'  shock  front 

(7) 

ii'y  a theory  developed  by  Yurkovich  (fluid  shock  code)  . The 
loading  pressures  for  the  bare  panel  problem  were  calculated 
v/itl.  the  fluid  drag  code  and  fluid  shock  code  and  input  into 
the  BR-1.^  (HR)  coda.  The  finite  element  code  does  not  have 
the  capability  of  calculati'''n  the  response  of  a wall  backed 
by  foam. 

The  second  approach  was  based  on  a two-dimensional  finite 
difference  coda  (AFTON)  together  with  a finite  element  code 
(WONSAP) . TFTON  finite  difference  calculations  were  used  to 
treat  the  bar.,  panel  cuse  as  well  as  the  foam  backed  panel 
case. 

The  threats  and  impact  configurations  selected  for 
computer  analysis  are  outlined  in  Taolo  2. 

5.1  THE  FLUID  DRAG  CODE  (f’DC) 

5.1.1  scription  O'f  the  Code 

The  fluid  draq  code  (FDC)  is  described  in  Reference 
3.  I c provides  a technique  for  calculating  the  fluid  pressure 
and  velocity  as  funorions  of  time  at  specified  points  within  a 

G9 


TABLE  2.  CALCULATION  MATRIX 


Panel  Considered 


Shot  No. 

Sphere 

Size 

(mm) 

Impact 

Velocity 

(km/s) 

Wall 

Material 

(Al) 

Foam 

Thickness 

(mm) 

BRl-A 

AFTON/ 

NONSAP 

FT5 

11.1 

1.53 

2024-T3 

0 

Entrance 

Entrance 

FT26A 

11.1 

1.52 

2024-T3 

38 

- 

Entrance 

FT2 

11.1 

1.73 

7075-T6 

0 

Entrance 

- 

FT  6 

14.3 

1.52 

2024-T3 

0 

Entrance 

Entrance 

FT8 

14.3 

1.57 

2024-T3 

38 

- 

Entrance 

FT33 

14.3 

1.83 

7075-T6 

0 

Side 

- 

tank  penetrated  by  a tumbling  projectile.  The  code  is  very 
efficient  and  uses  only  seconds  of  computer  time.  The  calcu- 
lated pressures  are  sufficiently  accurate  for  most  purposes 
as  long  as  certain  assumptions  made  in  the  code  are  not 
violated.  These  assumptions  are: 

1.  The  projectile  is  subsonic. 

2.  The  cavity  does  not  influence  the  pressure 
field . 

3.  Boundaries  of  the  fluid  are  either  perfectly 
rigid,  free,  or  perfectly  transmitting. 

4.  The  fluid  \/  1 \me  is  a rectangular  solid. 

5.  The  px'o.-iectile  drag  coefficient  is  inde- 
pendent of  projectile  velocity  (provision 
is  made  for  changing  drag  coefficient  due 
to  projectile  tumbling  and  case  stripping, 
but  these  were  not  used  for  this  work). 

In  previous  work^^^  this  code  was  used  to  calculate  the 


pressure  induced  by  tumbling  armor  piercing  (AP)  projectiles. 
The  side  panels  of  the  tank  were  so  far  from  the  trajectory  that 
their  influence  could  usually  bo  ignored.  Bullet  tumbling  and 
jacket  stripping  in  the  fluid  could  not  be  measured  very 
accurately;  however,  transducers  recorded  fluid  pressure  at 
typically  five  positions  in  the  tank.  Using  the  pressure  data, 
the  nine  parameters  in  the  FPC  which  pertain  to  bullet 
tumbling  and  stripping  wore  varied  to  bring  the  computed 
pressures  into  agreement  with  the  measured  pressures.  These 
parameter  values  could  then  be  used  to  compute  pressure  else- 
where in  the  target  tank. 

In  the  present  application,  no  tuiiibling  or  stripping  of 
the  projectile  occurs,  and  the  values  of  all  input  parameters 
were  uniquely  determined  from  measurements.  The  value  of  the 
drag  coefficient  (Cq)  and  the  initial  projectile  velocity  (u^) 
were  determined  from  the  experimental  data,  as  discussed  in 
Section  4.4, 

The  coordinate  system  defined  for  the  FDC  consisted 
of  an  origin  at  the  lower  right  hajid  corner  of  the 
entrance  panel;  + x pointed  across  the  tank,  + y was  up,  and 
+ 2 was  into  the  tank. 

The  volume  occupied  by  fluid  was  O^X'-x  , 0<y<y  , and 

0<z<2  , where  x , y , and  z were  C.850,  0.798,  and  0.980  m, 
c c c o 

respectively. 

Two  minor  modifications  were  made  to  the  FDC.  First, 
the  IMAGES  subroutine  was  modified  to  allow  boundary  conditions 
pertinent  to  the  UDBI  test  tank.  Pressure  reflections  from 
tank  boundaries  are  accounted  for  by  up  to  27  images  of  the 
projectile  moving  in  an  infinite  fluid.  In  the  original 
version  of  the  code,  all  images  were  of  the  same  type, 
depending  on  whether  the  boundaries  were  rigid  (+1  image) , 
free  (-1  image)  or  transmitting  (0  image) . In  the  present 
version,  the  values  of  the  ASIGN  matrix  are  altered  to  describe 
a tank  with  rigid  boundaries  at  x = 0,  and  at  y = 0..  free 
surfaces  at  x = x^,  y = y^,  and  z = z , and  a transmitting 
boundary  at  z = 0. 
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The  other  modification  necessary  to  the  code  was  a 
provision  to  calculate  and  print  the  input  data  in  the  format 
required  for  the  BR-lA  (HR)  code.  The  required  input  data  are 
just  twice  the  pressures  as  calculated  by  the  FDC  on  the 
z =0  plane.  It  was  determined  that  the  code  could  not 
calculate  pressure  at  z = 0,  so  the  entrance  panel  pressures 
were  calculated  for  z = 0.5  mm.  The  code  modifications  were 
inserted  into  the  MIRROR  subroutine  to  calculate  pressure  and 
into  the  HRAM  subroutine  to  write  data.  The  modification  made 
to  the  FDC,  as  well  as  sample  input  data,  are  listed  in 
Appendix  B. 

5.1.2  Validation  of  Pressure  rredictcd  by  the  FDC 

Validation  of  the  FDC  was  accomplished  by  comparing 
pressures  predicted  by  the  code  with  pressures  measured  by 
transducers  or  predicted  by  the  AFTON  code.  In  assessing  the 
code,  no  input  parameters  were  altered  to  improve  agreement 
with  the  measured  pressure  profiles. 

The  best  pressure  transducer  data  were  obtained  from 
stations  in  the  rear  half  of  the  target  tank.  For  Shot  FT33, 
the  computed  and  measured  pressure  pulses  for  the  T4  transducer 
are  compared  in  Figure  53.  The  transducer  was  situated  740  mm 
from  the  entrance  panel  and  100  mm  above  trajectory.  The  os- 
cilloscope record  was  filtered  with  a 30  kHz  (30  ys  rise  time) 
cut-off  filter.  Thus,  the  pulse  rise  time  may  have  been  lengthened 
and  reduced  slightly.  The  record  shown  in  the  figure  has  been 
smoothed  by  hand  to  eliminate  pulses  shorter  than  '^20  ys  duration. 
Negative  values  of  FDC  predicted  pressure  are  plotted  as  zero. 

The  agreement  of  the  two  pressure  profiles  is  very 
satisfactory.  The  arrival  time  of  the  predicted  pulse  agreed 
with  that  of  the  measured  pulse  better  than  the  estimated  measure- 
ment error  (±10  ys)  . Both  pulses  indicate  that  the  peak  pre.ssure 
is  attained  within  ^^>25  ys  after  first  arrival,  and  after  the  peak 
the  pressure  decays  rapidly.  After  about  50  ys  the  decay  be- 
comes less  rapid.  Over  two- thirds  of  the  total  impulse  is  con- 
tained in  the  more  slowly  decaying  regime,  which  lasts  for  about 
150  ys.  The  impulse  delivered  by  the  two  pulses  agree  within  the 
measurement  accuracy  (±15%). 
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Figure  53.  Comparison  of  pressure  computed  by  FDC  and  pressure 
measured  by  a transducer  in  Shot  FT33  at  a point 
located  at  100  mm  above  trajectory  740  mm  from  the 
entrance  panel. 

A spike  which  occurs  at  700  ws  in  the  FDC  calculation 
is  nearly  absent  in  the  transducer  record.  Inspection  of  the 
calculated  fluid  velocity  components  reveals  that  the  spike  is 
due  to  a reflection  from  the  bottom  of  the  tank,  which  was 
modeled  as  a rigid  boundary.  In  the  FDC,  the  reflection 
propagated  through  the  cavity  to  arrive  at  the  transducer 
station.  In  the  actual  experiment,  of  course,  the  cavity  does 
not  propagate  waves.  Fortunately,  late  time  arrivals  such  as 
these  carry  little  impulse  and  are  not  of  great  concern  in 
assessing  uhu  usefulness  of  the  FDC.  The  effect  of  similar 
cavity-related  errors  in  the  FDC  treatment  will  be  even  smaller 
at  the  panels. 

A similar  comparison  of  pressure  profiles  for  the  T3 
transducer  in  Shot  FT33  (at  z = 430  mm,  100  mm  above  trajectory) 
is  shown  in  Figure  54 . Consideration  of  the  indicated  zero 
shift  in  the  transducer  data  brings  the  two  profiles  into  good 
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Figure  54.  Same  as  Figure  53,  except  for  a point  0.43  m 
behind  the  entrance  panel. 

agreement.  The  same  general  features  noted  previously  apply 
here,  too.  On  the  computed  profile,  points  A,  B,  and  C mark, 
respectively,  the  arrival  of  reflections  from  the  bottom  surface, 
the  free  (top)  surface,  and  the  side  surfaces.  (Contrary  to  the 
appearance  of  the  figure,  it  can  be  seen  from  examination  of  the 
printouts  that  the  arrival  from  the  sides  had  no  effect  on 
pressure  and  dP/dt  did  not  change) . 

Figure  55  illustrates  the  agreement  between  the  FDC 
and  pressure  data  for  the  transducer  placed  nearest  the  entrance 
panel.  The  data  from  the  Tl  station  (z  = 133  mm,  220  mm  above 
trajectory)  in  Shot  FT5  are  used.  Again,  arrival  times  agree 
extremely  well.  When  the  transducer  record  comes  back  on  scale, 
it  is  in  satisfactory  agreement  with  the  FDC  pressure.  The 
agreement  is  improved  if  the  measured  spikes  at  230  and  290  ms 
are  regarded  as  instrument  noise.  This  is  probably  the  case, 
since  they  appear  on  almost  all  Tl  transducer  traces.  The 
accuracy  of  the  FDC  near  the  entrance  panel  can  be  assessed  by 
comparison  with  the  results  of  the  AFTON  calculations. 
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Figure  55.  Same  as  Figure  52,  but  for  a point  220  mni  above 
trajectory  and  133  mm  behind  entrance  panel. 

Figures  56  and  57  document  two  such  comparisons  within 
100  mm  of  the  panel.  The  FDC  does  not  model  shock  wave 
formation,  and  hence  fails  to  reproduce  the  leading,  peak  on  the 
AFTON  profiles.  The  AFTON  pressure  is  overall  slightly  higher 
because  in  the  AFTON  model  of  the  FT5  event  the  projectile 
enters  the  fluid  at  1.55  km/s,  whereas  in  the  FDC  model,  it 
enters  at  1.24  km/s.  The  general  agreement  between  these  two 
totally  different  numerical  treatments  implies  that  both  suc- 
cessfully model  the  effects  of  drag  forces  in  the  fluid. 

5.1.3  Results  of  Calculations  With  FDC 

The  FDC  was  very  helpful  for  obtaining  a qualitative 
understanding  of  several  phenomena  in  the  target  tanks.  For 
example,  most  transducer  records  showed  an  abrupt  drop  to  zero 
pressure  about  200  ys  after  first  arrival.  This  was  shown  to 
be  due  to  the  reflection  of  a pressure  wave  from  the  free 
surface.  It  was  also  shown  that  in  general  P>>PqU^  in  the  flow, 
so  that  the  orientation  of  a transducer  will  not  effect  its 
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reading.  This  was  confirmed  experimentally. 

Lundstrom  suggested'  ' that  negative  values  of  cal- 
culated pressure  be  set  equal  to  zero  for  computational 
purposes.  This  was  done  for  the  input  to  the  BR-lA  (HR)  calcu- 
lations. Lundstrom  has  suggested that,  at  least  near  the 
trajectory,  negative  values  of  pressure  may  signify  cavitation. 
In  Figure  58,  the  observed  cavity  at  z = 133  mm  is  compared 
with  regions  of  (y,  t)  space  in  which  the  Lundstrom  code 
predicts  negative  pressures. 
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Figure  58.  Comparison  of  experimental  and  predicted  cavitation 
at  z = 133  mm  for  Shot  FT5.  (Horizontal  lines 
show  region  in  which  calculated  pressure  is  less 
than  zero) . 
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The  initial  agreement  with  the  projectile  cavity  is  rather  good. 
However,  the  code  also  predicts  a cavitation  wave  from  the  upper 
free  surface  which  was  not  observed  within  the  space  visible  to 
the  trajectory-viewing  framing  camera. 

Figures  59  and  60  show  examples  at  front  panel  pres- 
sures calculated  for  FT5  and  FT5B.  It  can  be  seen  that  in  the 

central  regions  of  the  tank,  peak  pressures  are  approximately 
2 

7 MN/m  (1000  psi)  but  drop  by  a factor  of  10  in  about  300  ps. 
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Figure  59, 


Pressure-time  history  predicted  at  selected  points 
on  trajectory  for  impact  conditions  of  Shot  FT5. 
The  numbers  identify  pressure-history  points.  The 
coordinates  of  the  pressure-history  points  can  be 
found  from  Figure  64. 
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Figure  60.  Pressure-time  history  for  selected  points  above  and 
below  trajectory  for  impact  conditions  of  Shot  FT5. 
The  coordinates  of  the  pressure  history  points  are 
given  in  Figure  64. 

5.2  THE  FLUID  SHOCK  CODE 

The  shock  pressure  in  a hydrodynamic  ram  event  loads  a 
small  portion  of  the  entrance  panel  impulsively  during  the  first 
100  us  after  impact.  Pressure  profiles  behind  a shock  front  are 
extremely  difficult  to  measure,  but  using  the  theory  described 
by  Yurkovich^ , they  can  be  calculated  from  measurements  of 
the  shock  front  velocity.  His  treatment  uses  a nondimensional 
form  for  the  continuity  equations  and  the  equations  of  state, 
and  is  based  upon  the  work  of  Bach  and  Lee'  ' . For  this  study, 
only  the  portion  of  Yurkovich's  work  which  treated  the  pressure 
behind  a shock  front  is  described. 

The  conservation  equations  and  the  equation  of  state  are 
written  in  nondimensional  form.  The  following  definitions 
are  needed; 
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Shock  Mach  Number;  M = 


(t) 

c 


Position; 


r = 


®o 


Velocity; 


^ u Tty 

^ ^o 


Density; 


p(J,„  ) . £(^ 


Pressure; 


5,B,„  , = 

° ' (t) 

o 


(24) 


In  the  above  equations,  subscript  s^,  indicates  quanti- 
ties along  the  axis  of  symmetry,  c is  the  sound  velocity 
(1495  m/s  in  water)  and  a bar  designates  a non-dimensional 
quantity.  Using  these  definitions,  a profile  of  the  pressure 
behind  the  shock  front  can  be  described.  The  equation  for  the 
pressure  requires  the  non-dimensional  density  immediately  behind 
the  shock  front,  To  evaluate  this  quantity,  a non-dimensional 

form  of  the  adiabatic  equation  of  state  is  combined  with  the 
Hugoniot  relations  to  give; 


1 

n 


(25) 


when  n for  water  is  7.0  and  total  symmetry  was  assumed.  The 
general  equation  for  the  pressure  at  any  point  behind  the  shock 
is  then  given  by 
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c — , 


(q+2) 


1 - / - - \ 
^2  P3  I 1+  r [ {q+2)  (^nr)  -l]j* 


(26) 


_ q+2^ 

+ ^ P.  [ 2-  r J [ (q+2)  (inr)  -1]  ^ +1 

(q+2)^  ^ \ I 


The  non-diraensional  pressure  terms  for  Equation  (25)i.-are  given 
by; 


P,  = 1 - 

Pi 


(27) 


P2  = (5^-1)  (1-H)  + e (G^  + Mg  1^), 

o s_ 


(28) 


P,  = 


U^  K (l+Uj^(H-2)  - 0]  - 0 Mg  (Uj^H)  , 

° ®o 


(29) 


-2  2 


(30) 


where ; 


^ ^8 


S 


(31) 


and 


BM 


H 


s 


3M_  Pi 


(32) 


For  the  case  of  a strong  shock,  the  general  form  of  the  equation 


can  be  simplified  to; 


P = ?!  + (u^  + 8-1).  (33) 

Yurkovich  made  the  assumption  that  the  shock  would  be  strong . 
The  simplified  initial  conditions  led  to  the  equation 

(34) 


+ (1-r 


BM, 


(Pj^)  (q+2) 


DM. 


The  calculation  of  the  pressures  from  the  above  equations  is 

dependent  upon  the  measurement  of  the  shock  front  Mach  number 

M , and  the  determination  of  the  parameter  8. 
s 

For  this  work,  8 was  found  by  using  Equation  (31)  and 
assuming  that  the  shock  front  position  as  a function  of  time  is 
described  by 


r.  « kt'^  (35) 

w 

Yurkovich  selected  the  form  of  Equation  (34 ) for  convenience  and 
because  it  described  the  shock  waves  produced  by  hypervelocity 
impact  in  previous  experiments.  Yurkovich  used  ij;  » 0.8  as 
descriptive  of  water.  The  data  for  Shot  FTA4  are  fit  by 
Equation  (35)  with  the  parameters  k » 2.54  mm  and  (j;  » 0.901 
(for  t in  us)  or 

r^  » 2.54t°*^°^  (36) 

s 

A computer  routine  to  perform  the  shock  pressure  calcu- 
lations was  written  and  is  listed  in  Appendix  B.  The  results 
of  three  methods  of  calculating  the  pressure  are  illustrated 
in  Figure  The  non-dimensional  pressures  are  plotted  as  a 
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Figure  61.  Typical  pressure  profile  calculated  by  three  methods. 

function  of  the  non-dimensional  radius  of  the  shock  front.  For 
the  illustrated  pressures,  the  shock  front  Mach  number  was 
1.25  and  the  parameter  6 was  -0.10988.  From  Equation  (35)  , this 
corresponds  to  a shock  radius  of  16  mm  at  a time  of  7.74  us  after 
impact . 

The  three  pressure  solutions  converged  at  the  shock  front 
2 

to  384  MN/m  (55.7  ksi) . Near  r = 0 the  computed  pressure 
significantly  different  for  the  three  methods.  However,  it 
should  be  noted  that  the  entrance  hole  extends  from  r = 0 to 
0.347.  No  pressure  can  be  exerted  on  a wall  over  this  region. 

Ni.  ne  of  the  three  methods  for  calculating  the  pressures  on  the 
wall  near  the  entrance  hole  takes  into  account  the  penetrating 
projectile. 

The  Yurkovich  pressure  equation  (Equation  (34))  was  used 
to  generate  the  impulse  delivered  to^  an  element  of  the  entrance 
panel  imiaediately  after  penetration.  The  scheme  for  calculating 
the  impulse  is  shown  in  Figure  62.  The  force  on  each  region  is 
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figure  62-  Scheme  used  in  calculating  the  impulse  delivered  by 
the  shock  pressures  to  the  entrance  panel.  No 
forces  are  applied  to  the  entrance  hole  (shaded 
area) . Forces  on  the  remainder  of  the  region  are 
calculated  from  pressure  averages. 


computed  as  the  average  pressure  (P^  + ^i-1^^^  times  the  area  of 


the  region  A^.  The  area  is  given  by 

A . = li  ( R . ^ - R . ^ , ) 

1 1 1 -1 

th 

and  the  force  on  the  i^  region  is  given  by 

2- 

P . = p U P . A. 

1 ^O  S 1 1 


(37) 


(38) 


The  average  pressure  for  the  i section  is  read  from  pressure 
profile  plots  similar  to  those  shown  in  Figure  61.  Total  forces 
were  calculated  at  several  different  times  after  impact.  The 
forces  after  impact  for  Shots  FTA4  and  FT5F<  are  shown  in 
Figure  63.  The  forces  were  integrated  to  obtain  impulse  values 
which  were  input  into  the  BR-lA  (HR)  structure  code- 
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Fiyuro  (>3,  Plot  of  f-'rco  vs,  tiuio  for  pressure  calculated  at 
entrance  panel  for  Shot  FTSB.  Impulse  was 
calculated  from  this  plot  numerically. 

5.3  BR-IA  (HU)  STRUCTURAL  RESPONSE  CODE 

Deflection  of  bare  entrance  and  side  panels  was  numori- 

caily  computed  using  the  BR-lA  (HR)  structural  response  code. 

(15) 

The  basic  code,  BR-IA,  was  developed  by  Northrop  Corporation 
to  model  blast  loading  of  aircraft  structures.  For  application 
to  hydrodynamic  ram  problen^s,  the  code  was  modified  (HR)  at  the 
Naval  Post  Graduate  School 

Considerable  effort  was  expended  in  developing  a correct 
version  of  the  BR-IA  (HR)  code.  The  original  (HR)  modifications 
were  made  on  an  IDfi  360  version  of  the  BR-1  code.  After  these 
modifications  were  made  to  the  BR-lA  code  and  a sample  problem 
executed  on  the  CDC  6600,  the  results  did  not  agree  with  those 
of  Reference  16.  An  error  in  the  CDC  6600  BR-lA  code  was 
finally  uncovered.  The  statement  labeled  3431  in  the  XPORCE 
routine  should  read  3431  JJ  ■ JJ  + 2. 

The  HR  modification  to  the  code  were  based  on  one- 
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dimensional  piston  theory  which  calculates  the  pressure  on  a 
piston  as 


P = + pc  (Uj  - w)  , (39) 

where  p is  the  net  pressure  on  a boundary  wall,  p.  is  the  inci- 
dent fluid  pressure  (e.g.,  the  pressure  as  if  the  wall  did  not 
affect  the  fluid) , u is  the  incident  fluid  velocity  normal  to 
the  wall,  and  w is  the  wall  velocity. 

The  assumption  that  P^  = peu  leads  to 


2P^  - pew. 

(40) 

• 

pc  (2Uj-w) . 

(41) 

The  physical  basis  of  these  forms  of  piston  theory  is 
discussed  in  Reference  22.  All  are  exact  for  planar  incident 
waves  at  normal  incidence.  Equation  (40)  is  exact  for  an 
arbitrary  incident  wave  on  a plane  rigid  wall,  and  Equation  (41) 
is  exact  for  an  arbitrary  incident  wave  on  a plane  free  surface. 

For  nearly-free  surface  boundary  conditions  Equation  (40 
should  give  better  solutions  than  Equation  (39) . Thin  aluminum 
walls  are  expected  to  behave  more  like  free  surfaces  at  early 
times  . nd  like  rigid  boundaries  as  membrane  strains  develop.  On 
the  other  hand,  in  the  motions  which  are  considered  in  this 
report,  it  turns  out  that  P^>>pcu^.  Thus,  the  use  of  Equation 
(40)  will  lead  to  greater  driving  stresses  than  the  use  of 
Equation  (41) . Since,  in  general,  use  of  the  BR-lA  (HR)  code 
leads  to  underprediction  of  displacement,  we  have  used  the  form 
of  piston  theory  of  Equation  (40)  throughout.  This  same  proce- 
dure was  used  in  Reference  9 for  the  same  reason.  With  this 
form  of  piston  theory,  the  input  required  for  the  BR-lA  (HR)  is 
twice  the  incident  pressure  P^,  which  would  be  present  if  there 
were  no  fluid-wall  coupling.  This  pressure  was  provided  by 
the  FDC  and  fluid  shock  programs  described  in  the  preceding 
sections. 


5.3.1 


Calculations  of  Entrance  Panel  Motion 
The  elements  used  to  model  one-half  of  the  front  face 
of  the  tank  are  shown  in  Figure  64.  circled  numbers  designate 
the  plate  elements  and  the  uncircled  numbers  designate  the 
joints.  The  projectile  impact  occurs  at  joint  32.  The  line 
between  joints  2S  and  36  is  a line  of  synm\etry.  The  bottom 
and  top  edges  of  the  tank  are  assumed  simply  supported;  how- 
ever, the  slant  side  is  not  as  easily  modeled.  The  allowable 
constraints  provided  by  the  code  are  in  the  global  coordinate 
system,  and  it  is  desirable  to  let  the  slant  edge  rotate  about 
a line  parallel  to  the  edge.  This  was  done  by  allowing  side 
joints  to  rotate  about  both  the  x and  y axes. 


Figure  64  . Finite  element  model  of  left  side  of  a fuel  tank 
entrance  panel.  Impact  occurs  at  joint  32.  The 
right  side  of  the  tank  wall  is  treated  as  synm\etric. 
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One  very  important  shortcoming  of  the  program  is  its 
inability  to  handle  the  initial  static  deflection  due  to  the 
hydrodynamic  pressure  in  the  tank.  It  was  initially  thought 
that  this  could  be  done  by  applying  the  hydrostatic  pressure 
as  a step  function,  running  the  problem  to  equilibrium  and  then 
appl^iing  the  shock  and  drag  pressures.  However,  due  to  the 
small  time  step  required  by  the  program  for  numerical  stability, 
this  was  not  possible  without  an  unacceptible  use  of  computer 
time.  An  initial  attempt  at  this  problem  using  100  time  steps 
(total  time  of  0.0014  seconds)  required  570  seconds  of  CP  time 
and  over  an  hour  of  10  time  on  the  CDC  6600.  A run  to  equili- 
brium might  take  ten  times  this  long.  The  efficient  solution 
of  this  problem  would  require  modification  of  the  program  to 
provide  a solution  to  the  initial  static  problem. 

An  additional  weakness  in  the  program  is  the 
necessity  of  providing  input  pressures  which  are  constant  over 
an  element.  This  makes  modeling  of  the  initial  impact  and 
early  drag  forces  rather  arbitrary  since  the  pressure  gradients 
are  very  large.  This  problem  could  be  ameliorated  but  not 
eliminated  by  using  smaller  elements. 

5 . 3 . 1 . 1 Case  of  11.1  mm  Diameter  Projectile  Striking 

2024-T3  Entrance  Panel  at  1 ,'5  km/i" 

This  case  corresponded  to  experiments,  FT5,  FT5A, 
FT5B,  and  FTA4  and  was  also  calculated  with  the  AFTON/NONSAP 
codes  (see  Appendix  A) . Samples  of  the  pressures  calculated 
with  the  FDC  and  input  into  the  BRl-A  (HR)  code  are  given  in 
Figures  65  and  66.  The  running  time  was  typical,  431  s CP  and 
1993  s 10. 

It  is  apparent  that  the  discrepancy  between  pre- 
dicted and  observed  panel  displacement  is  very  severe.  In  the 
calculation,  the  peak  displacement  is  reached  in  about  600  ys, 
and  is  15  percent  of  the  eventual  peak  displacement  observed 
for  that  experiment.  The  calculation  also  fails  to  show  the 
relatively  sharp  increase  in  displacement  which  propagates 
away  from  the  impact  point.  Thus,  the  computed  panel  dis- 
placement profile  is  neither  in  qualitative  nor  quantitative 
agreement  with  observation. 


Figure  65.  Comparison  of  panel  profiles  predicted  by  BR-lA  (HR) 
code  and  measured  during  Shot  FT5B.  Impact  condition 
is  11.1  mm  diameter  sphere  striking  1.5  mm  thick 
2024-T3  aluminum  panel  at  1.5  km/s.  BR~1A  (HR) 
was  driven  by  pressures  calculated  with  FDC  and 
shock  code. 


Figure  66 . Out-of-plane  displacement  at  nodes  of  calculation 
shown  in  Figure  65  , at  0.89  ms  after  impact. 


The  outward  velocity  of  the  panel  at  node  16,  which 
is  202  nmi  away  from  the  impact  point,  is  shown  in  Figure  67  , 

The  panel  begins  to  move  out  at  80  us  after  impact;  this  motion 
is  due  to  the  affects  of  the  early-time  shock,  which  have  pro- 
pagated through  the  aluminum.  The  pressure  wave  in  the  fluid 
arrives  at  about  130  us.  At  that  time  the  panel  accelerates 
outward.  At  times  later  than  about  200  us,  the  panel  decelerates, 
and  by  600  us  the  panel  has  stopped  moving.  This  behavior  is 
completely  contrary  to  observation. 

The  deceleration  of  the  panel  can  only  be  due  to  an 
inward-directed  force.  The  strains  in  the  panel  at  these  times 
are  too  small  to  provide  such  a force.  The  restraining  forces 
arise  from  the  piston  theory  embodied  in  the  DR-lA  (HR)  code. 

The  piston  theory  term  for  this  node,  2P-ticw,  is 
qraphed  in  Figure  68.  (The  incident  pressure  can  be  seen  in 
Figure  59)  . Note  chat  this  term  is  generally  negative.  This 
is  because  the  incident  pressure  falls  off  rapidly  from  its 
initial  maximum  value.  Physically,  such  an  incident  wave  would 
cause  the  front  wall  to  spall  away  from  the  fluid.  However, 
the  piston  theory  does  not  allow  cavitation,  and  so  a novjative 
pressure  is  exerted  on  the  panel  causing  it  to  decelerate  and 
stop.  Previous  investigators  have  also  noted  this  phenomenon, 
both  in  calculations and  experiments ^ . 

Another  cause  for  error  in  the  BR-lA  (HR)  treatment 
of  entrance  panel  motion  lies  in  the  derivation  of  the  incident 
pressures  themselves.  It  is  assumed  that  the  pressure  incident 
on  any  element  of  the  panel  is  independent  of  the  motion  of 
other  elements  of  the  panel.  This  is  clearly  not  the  case. 

Motion  of  the  panel  in  the  region  near  the  entrance  hole  will 
attenuate  fluid  pressure  waves  before  they  reach  more  distant 
elements.  This  is  particularly  evident  in  the  case  of  the  shock 
wave.  As  discussed  in  Section  IV,  the  observed  shock  fronts 
are  not  hemispherical,  which  indicates  that  the  peak  pressure 
at  the  front  decays  more  rapidly  along  the  panel  than  in  tlie 
fluid;  this  behavior  is  not  mimicked  by  the  calculations. 
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07.  Panel  velocity  at  node  16  (see  Figure  64)  for 
calculation  shown  in  Figure  65. 


Figure  68.  Piston  theory  pressure  at  node  16  for  the  calcula 
tion  shown  in  Figure  65. 


Note  that  those  two  sources  of  error  are  in  opposite 
directiox^s.  The  failure  to  treat  cavitation  leads  to  reduced 
displacement,  while  the  failure  to  treat  incidexxt  wave  atteixua- 
tion  leads  to  increased  displacement. 

In  order  to  evaluate  the  BR-1  code  without  the 
uncertainties  iixtroduced  by  piston  theory,  a calculation  of  the 
FT5  experiment  was  carried  out  usinvj  pressure  data  yenorated  by 
the  AFTON  code.  The  pressure  was  calculated  at  points  1 nun 
behind  the  entrance  panel,  which  was  modeled  by  mass-loads 
only  (e.>;].,  no  tensile  stren<jth)  . The  pressure  data  supplied 
by  CRT  were  iixtosjraced  and  ixxterpolated  by  haiid  to  derive  input 
data  for  a BK-IA  (run  without  the  HR  option).  The  results 
are  shown  in  Fiqure  i'"'*. 


Figure  69 . Comparison  of  BR-IA  calculation  with  experimental 
data  for  case  of  11.1  nun  diameter  sphere  striking 
1.5  mm  thick  2024-T3  panel  at  1.5  kiVs.  BR-IA  was 
driven  by  pressures  calculated  by  AFTON  code. 
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The  displacements  predicted  by  this  calculation  are 
up  to  four  times  yreater  than  those  from  the  previous  oaloulat  ioii . 
There  has  also  been  a significant  qualitative  improvement  in 
the  profile.  At  940  ps  the  panel  is  continuing  to  move  outward, 
and  it  appears  the  final  displacements  will  be  about  twice 
the  values  at  that  time.  The  regions  of  relatively  sharp 
increase  in  displacement  at  about  250  mm  and  90  mm  are  qualita- 
tively reproduced.  Quantitative  accuracy  of  the  calculated 
displacement  at  900  us  is  still  poor;  the  peak  value  is  low  by 
a factor  four.  It  appears  that  beyond  the  inunediate  entrance 
region,  final  displacement  will  be  about  50'i  too  low. 

The  displacements  are  .tlso  considerably  less  than 
those  of  the  corresponding  NON?AP  calculation  (Figure  .V^l)  , At 
too  us  the  discrepancy  is  rinly  about  15‘f. , but  in  the  onti  ance 
region,  it  is  a factor  3.  At  I ms  the  agreement  is  improvt.Hi 
because  the  NONSAP  panel  has  begun  to  rebound,  and  the  AR-IA 
panel  has  not . 

The  difference  between  these  two  calculations  is 
apparently  duo  to  the  number  of  elements.  The  NONSAP  calcula- 
tion used  56  annular  elements.  The  BR-lA  used  32.  However,  a 
horizontal  radius  in  the  UR-IA  calculation  only  passes  through 
tour  elements,  the  first  of  which  is  76  mni  wide.  Otven  the 
coarse  zoning,  the  degree  of  disagreement  between  these  caleu- 
lations  is  not  unreasonable.  It  is  also  clear  that  the  HR-lA 
cannot  treat  the  region  immediately  adjacent  to  impact  witfiout 
a great  increase  in  zone  density,  and  hence  running  time. 

One  run  was  carried  out  in  which  the  usual  driving 
pressures  were  inserted  into  the  BR-IA  (HR)  code,  but  the  (HR) 
option  was  not  exercised.  The  particular  input  pressures  used 
were  those  from  Shot  FT2,  which  was  at  a 12%  higher  impact 
velocity  than  the  shot  considered  above.  As  expected,  neglect 
of  the  piston  theory  yielded  dramatic  displacements.  Zones 
around  the  entrance  panel  failed  at  once.  This  calculation  was 
run  out  to  145  us.  At  that  time  the  peak  deflection  was  over 
75  mm.  Zones  adjacent  to  the  impact  site  were  moving  out  at 
250  m/s.  Thus,  neglect  of  piston  theory  is  not  a viable  means 
of  bringing  BR-IA  calculations  into  agreement  with  observation. 
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5 . 3 . 1 . 2 Case  of  11.1  nun  Diameter  Pro:]ectilo  Striking  7075-T6 

Entrance  Panel  at  1.7  3 kni/s 

This  case  corresponds  to  Shots  FT2  and  FT3.  In 
both  shots,  the  entrance  panels  failed  catastrophically. 

The  FDC  calculation  was  carried  out  with  Cj^  * 0.356 
and  an  entrance  velocity  of  1.26  km/s.  The  shock  loading  was 
obtained  by  scaling  up  the  impulse  loads  from  the  previous 
calculation  by  the  ratio  of  impact  velocities. 

No  failure  in  any  zone  was  predicted  by  the  code. 

The  computed  wall  deflections  were  once  again  minimal.  A hori- 
zontal profile  across  trajectory  is  shown  in  Figure  70 . 


Figure  70.  Displacement  profile  predicted  for  case  of  11.1  nmi 
diameter  sphere  striking  1.5  mm  thick  7075-T6  panel 
at  1.73  km/s.  Contrary  to  experiment,  no  failure 
was  indicated  by  calculation. 
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5. 3.1. 3 


Case  of  14.3  itmi  Diameter  Sphere  Striking  2024-T3 
Panel  at  1.52  km/s. 


This  case  corresponds  to  Shots  FT6  and  FTA9 . A 
problem  was  encountered  in  the  shock  analysis  for  this  run.  The 
decay  parameter,  n,  determined  along  trajectory,  was  found  to 
be  0.988.  This  was  apparently  caused  by  the  motion  of  the  large 
projectile  driving  the  shock.  The  shock  fronts  calculated  by 
the  fluid  shock  code  with  this  value  of  n decayed  slower  than 
1/r;  thus,  the  energy  deposited  in  the  entrance  panel  by  the 
shock  wave  increased  with  time.  This  is  a non-physical  effect 
caused  by  the  fact  that  attenuation  of  pressure  at  outlaying 
zones  by  close-in  zones  cannot  occur  in  the  calculation.  To 
keep  some  degree  of  realiti-  in  the  input,  the  shock  loading  was 
arbitrarily  terminated  at  50  ys.  The  calculaiad  wall  displace- 
ments, shown  in  Figure  70,  reveal  the  effect  of  the  intense 
shock  loading.  However,  no  failure  was  indicated  by  the  calcu- 
lation, even  though  massive  failure  occurred  in  the  actual 
experiment. 

5.3.2  Calculations  of  Side  Panel  Motion 

The  elements  used  for  calculations  of  side  panel  motion 
are  shown  in  Figure  71.  The  nomenclature  is  the  same  as  for 
Figure  64.  The  boundaries  are  taken  as  simply  supported.  The 
panel  properties  were  chosen  for  7075-T6  aluminum. 

5 . 3 . 2 . 1 Case  of  14.3  mm  Diameter  Sphere  Striking  at  1.83  km/s 

This  case  corresponds  to  Shots  number  FT32,  FT33  and 
FTA9 . The  FDC  code  was  run  with  a coefficient  of  drag  equal  to 
0.334.  The  entrance  velocity  of  the  sphere  in  the  fluid  was 
taken  as  1.56  km/s.  As  noted  in  Section  IV,  the  agreement  of 
predicted  and  experimentally-measured  fluid  pressure  in  the 
tank  interior  was  rather  good  for  this  case.  The  FDC  calcula- 
tion was  carried  out  to  1.25  ms,  at  which  time  the  panel  pressures 
were  negligible.  The  BR-lA  (HR)  calculation  was  carried  out  to 
2.5  ms . 

The  tank  displacements  at  trajectory  height  predicted 
by  the  calculation  are  shown  in  Figure  72.  Again  the  displace- 
ments are  far  less  than  chose  observed  for  Shot  FT33.  The 
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Figure  71.  Element  used  in  BR-IA  calculation  of  side  panel 
displacement. 
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Figure  72. 
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Comparison  of  side  panel  profile  computed  by 
BR-IA  (HR)  calculation  and  profile  measured  in 
experiment  FT33. 
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calculated  profiles  show  a shape  which  is  basically  correct, 
but  especially  at  later  times,  their  amplitudes  are  only  about 
15%  of  the  actual  values. 

The  failure  of  the  code  to  account  for  maximum  dis- 
placement is  again  apparently  due  to  its  neglect  of  cavitation. 
Figure  71  illustrates  the  computed  wall  velocity  at  a point  near 
the  center  of  the  panel.  The  loading  pressure  applied  to  this 
region  of  the  panel  rises  to  a peak  in  10  us  and  falls  to  half 
its  peak  value  in  30  us;  it  has  a second  arrival  at  '\^200  us 
after  first  arrival.  Thus,  it  turns  out  that  the  piston  theory 
term,  2P-pcw,  becomes  negative  within  100  us.  The  reflected 
arrivals  from  the  bottom  and  opposite  side  of  the  tank  are 
evident  in  the  velocity  history.  However,  due  to  the  retarding 
force,  by  1.1  ms  the  panel  has  no  outward  velocity.  In  con- 
trast, the  actual  panel  has  an  average  velocity  of  about  4 m/s 
between  0.9  and  1.9  ms  and  did  not  stop  until  'uO  ms. 

5.4  FINITE  DIFFERENCE  CALCULATIONS 

Another  computational  routine  used  in  this  hydrodynamic 
ram  study  was  the  AFTON  code  of  California  Research  and  Tech- 
nology (CRT) . It  is  a two-dimensional,  finite  element  wave 
propagation  code.  It  is  an  established  program,  chosen  for 
this  work  because  it  has  been  used  on  a large  number  of  impact 
and  penetration  problems  similar  to  hydrodynamic  ram. 

5.4.1  .50  Caliber  Projectile  Impacts 

As  part  of  the  present  program,  several  WAVE-L  (a  code 

similar  to  AFTON)  calculations  of  hydrodynamic  ram  events  were 

executed  . A study  of  .50  caliber  projectile  impacts  was 

(5)  . . 

completed  and  reported  separately  . A .50  caliber  projectile 
was  modeled  and  allowed  to  penetrate  an  axisymmetric  tank  through 
a flat  entrance  panel  and  exit  the  tank  through  a flat  exit 
panel.  The  cross-sectional  area  of  the  projectile  was  varied 
to  simulate  the  changing  coefficient  of  drag  observed  during 
tumbling.  Two  tank  models  were  used,  both  had  3.2  mm  aluminum 
entrance  and  exit  panels.  The  tanks  had  different  depths,  76 
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Figure  73.  Panel  velocity  at  node  27  (see  Figure  71^  of 
calculation  shown  in  Figure  72. 

and  381  nun.  The  76  nun  tank  model  was  used  for  code  verification 
and  for  an  investigation  of  the  effects  of  rigid  foam.  The 
381  mm  tank  model  was  used  for  the  principal  studies  where 
there  was  no  foam  in  the  tank.  An  illustration  of  the  pro- 
jectile penetrating  the  entrance  wall  is  shown  in  Figure  74 . 

As  the  projectile  penetrated  the  fluid,  the  code  calculated 
fluid  particle  velocity  and  pressures  and  tank  wall  motion. 

Fluid  cavitation  was  observed  as  well  as  wall  stresses  suffi- 
cient to  initiate  cracks. 

As  the  projectile  traversed  the  381  cm  tank,  it  was 
allowed  to  tumble  rapidly  and  was  fully  tumbled  when  it 
reached  the  exit  wall.  The  resulting  coefficient  of  drag  is 
shown  in  Figure  75  contrasted  to  that  postulated  for  a more 
slowly  tumbling  projectile.  For  comparison,  the  coefficient  of 
drag  of  a high  velocity  steel  sphere  in  water  is  also  shown. 
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Figure  74. 


Finite  element  model  of  .50  caliber  projectile 
perforating  the  entrance  wall  of  a fuel  tank. 


Rigid  foam  backing  to  the  exit  panel  was  also  modeled. 

The  equation  of  state  for  AVCO  Thermarest  rigid  foam  was  used. 

The  compression  of  the  foam  was  also  computed  with  the  finite 

difference  code.  A principal  characteristic  of  the  rigid  foam 

is  that  it  transfers  all  of  the  stress  that  it  receives  until 

that  stress  reaches  the  crushing  strength  of  the  foam  (approxi- 

2 

mately  0.17  MN/m  ).  The  foam  then  yields,  allowing  expansion 
of  the  fluid  and  subsequent  reduction  of  the  fluid  pressure. 
Studies  with  a one-dimensional  model  of  the  foam  indicated  that 
approximately  30  to  40  mm  of  foam  would  adequately  protect  a 
3.2  mm  thick  exit  panel.  From  these  studies,  it  was  recommended 
that  the  entrance  panel  be  thin  so  that  only  a small  amount  of 
energy  would  be  transferred  directly  to  the  panel  during  pene- 
tration. It  was  also  recommended  that  the  entrance  panel  have 
a large  surface  area  backed  by  foam  which  would  absorb  the 
energy  of  the  penetrating  projectile. 

5*4.2  Calculation  of  Hydrodynamic  Ram  Effects  of  High 
Velocity  Fragments 

For  the  high  velocity  fragment  calculations  a slightly 
different  calculational  technique  was  employed.  The  WAVE-L 
calculations  conducted  for  the  projectiles  presented  two  im- 
portant difficulties.  Firstly,  the  Lagrangian  formulation  of 
the  WAVE-L  code  did  not  model  the  high  velocity  flow  around 
the  projectile  very  well.  Secondly,  the  finite  difference 
formulation  was  not  an  efficient  way  to  model  the  wall  response. 

The  details  of  the  calculation  approach  and  principal 
results  are  given  in  Appendix  A. 

5. 4. 2.1  Outline  of  Approach 

Recognizing  that  numerical  solutions  of  impacts 
into  fuel  tanks  must  be  efficient  (inexpensive)  in  order  to  be 
useful,  a number  of  important  innovations  were  made  for  the 
high  velocity  fragment  calculations.  First,  a coupled  finite 
difference  (FD) /finite  element  (FE)  approach  was  implemented 
to  avoid  the  serious  problem  which  would  be  caused  by  tiny 
time  steps  if  the  plate  were  to  be  treated  in  a FD  code. 

Second,  a finite  difference  code  (AFTON)  which  has  characteris- 
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-tics  which  are  well-suited  to  the  fuel  tank  application  was 

adapted  and  demonstrated. 

Initially,  a 11.1  ram  diameter  sphere  at  1.55  km/s 

(Case  A,  corresponding  to  Shot  FT5)  was  calculated  with  AFTON 

until  the  sphere  penetrated  many  diameters  and  the  entrance 

2 

panel  pressures  decayed  to  the  0.1-1  MN/m  level  (a  time  of 
175  ys) . 

The  pressure-time  histories  along  the  entrance 
panel  provided  the  boundary  conditions  to  a FE  calculation 
(using  the  NONSAP  code)  of  the  panel,  which  was  then  carried 
out  to  3 ms. 

Comparisons  between  some  of  the  calculated  and 
measured  panel  displacements  for  this  configuration  are  shown 
in  Figure  76.  (Additional  comparisons  are  shown  in  Appendix  A) . 


HORIZONTAL  DISTANCE  FROM  IMPACT  PONT  (mm) 


Figure  76.  Comparison  of  AFTON-NONSAP  calculation  with  experi- 
mental data  for  case  of  11.1  mm  diameter  sphere 
striking  1.5  mm  thick  2024-T3  panel  at  1.5  km/s. 
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Figure  76  shows  that  at  early  times,  up  to  about  C.6  ms,  the 
computed  and  measured  profiles  are  in  approximate  agreement. 

In  fact,  it  can  be  shown  that  the  predicted  peak  displacement 
is  within  the  measurement  uncertainty  of  the  experimental  value. 
However,  apparently  due  to  the  premature  termination  of  driving 
pressures  in  the  calculation,  the  computed  profiles  rebound  at 
later  times,  when  the  actual  panels  continued  to  move  out. 

This  can  be  seen  in  the  figure  by  comparing  the  two  profiles 
for  1 ms  after  impact. 

Case  B (Shot  FT26) , in  which  the  panel  was  backed 
by  38  mm  of  foam,  was  then  run  to  175  us  on  AFTON.  The  entrance 
panel  pressures  calculated  in  Case  B were  orders  of  magnitude 
smaller  than  had  been  calculated  in  Case  A.  It  was,  therefore, 
clear  that  failure  would  not  occur  in  Base  B,  and  it  was  decided 
to  concentrate  on  an  evaluation  of  the  sensitivity  of  the 
NONSAP  finite  element  analysis  to  some  of  the  input  conditions 
and  assumptions. 

Consequently,  several  additional  NONSAP  calculations 
of  Case  A were  performed,  namely; 

(a)  A FE  solution  was  obtained  of  the  static  dis- 
tortion of  the  panel  under  a uniform  0.05  MN/m*  hydrostatic  load, 
representing  the  pre-impact  water  head  at  the  level  of  the  impact 
point  on  the  tank.  This  solution  produced  a maximum  displacement, 
at  the  center  of  the  panel,  of  5 mm.  This  was  about  half  the 
static  displacement  which  was  observed  in  the  test  panels.  This 
difference  indicates  that  the  test  panel  shape  (trapezoidal)  and 
edge  constraint  conditions  substantially  affect  its  static 
response. 

(b)  The  dynamic  FE  calculation  of  Case  A was 
repeated,  except  that  the  panel  was  initially  statically 
deformed  to  simulate  the  effect  of  the  static  hydrostatic  load 
on  the  dynamic  displacements.  The  minimal  differences  between 
this  repeat  and  the  original  Case  A analysis  indicated  that 
the  initial  static  distortions  in  this  case  were  not  signifi- 
cant to  dynamic  response. 
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(c)  Several  short  calculations  in  which  the 
inelastic  properties  of  the  aluminum  were  varied  did  not  produce 
changes  that  would  significantly  alter  the  calculated  results. 
These  trends  established  that  the  strength  modeling  in  NONSAP 
was  adequate. 

Based  on  the  experience  gained  from  the  preceding 
effort#  when  the  third  case  was  examined  (a  14.3  mm  diameter 
sphere  impacting  at  1.5  km/s)  two  considerations  led  to  a first 
attempt  of  running  a FE  calculation  (Case  C-1)  using  the  AFTON 
loading  obtained  from  Case  A,  scaled  by  the  ratio  of  the  sphere 
diameters.  First,  this  approach  would  be  very  attractive  for 
future  design  studies  of  impacts  onto  fuel  tank  walls  (with  or 
without  foam  backing) , since  a single  AFTON  FD  analysis  could 
be  used  as  input  to  several  plate  response,  FE  calculations, 
thereby  significantly  reducing  costs.  Second,  since  the  en- 
trance panel  thickness  doesn’t  scale  when  the  sphere  diameter 
is  changed,  this  procedure  slightly  overestimates  the  impulse 
delivered  to  the  front  panel.  Thus,  if  no  difference  between 
Case  C-1  and  Case  A was  observed,  the  ability  of  the  numerical 
technique  to  predict  front  panel  failure  would  be  suspect.  As 
it  turned  out,  the  resulting  stresses  and  strains  in  Case  C-1, 
as  stated  in  the  appendix,  were  sufficiently  different  from 
Case  A to  provide  convincing  evidence  that  the  code  solution 
could  differentiate  between  plate  response  to  various  loadings. 

The  scaled  impulses  used  in  the  foregoing  step  were 
expected  to  bo  overestimates  of  the  test,  since  the  actual  panel 
thicknesses  did  not  scale.  Therefore,  a PD  calculation.  Case 
C-2,  with  the  14.3  n\m  diameter  sphere  was  run  with  the  correct 
panel  thickness.  Peak  pressure  and  impulses  along  the  panel  in 
Case  C-2  were  found  to  be  only  slightly  smaller  than  the  scaled 
values  from  Case  A which  were  used  in  Case  C-1  (see  Figure  A-18 
in  appendix) . This  confirmed  the  general  validity  of  the  scaled 
loading  approach.  Comparisons  of  other  measured  data  (sphere 
displacement,  shock  velocity,  water  pressures)  were  excellent. 


103 


Original  plans  for  a second  foam-backed  panel 
calculation  were  abandoned  because  the  Case  B results  and  the 
experimental  data  showed  that  the  38  mm  thick  foam  was 
considerably  thicker  than  needed  to  prevent  fractures.  Thus, 
another  calculation  with  such  a thick  foam  layer  would  not  be 
expected  to  produce  any  additional  information.  The  remaining 
effort  was,  therefore,  expended  on  examining  refinements  of 
the  numerical  techniques  which  could  substantially  reduce  the 
cost  per  case  and  improve  accuracy  of  these  solutions . 

5. 4. 2. 2 Conclusions 

The  principal  conclusions  from  the  AFTON-NONSAP 
calculations  are  as  follows: 

1.  The  AFTON  code  is  capable  of  accurately 
generating  the  pressure  field  within  the  fluid  for  the  length 
of  time  that  loading  of  the  entrance  panel  is  significant. 

2.  AFTON-driven  NONSAP  calculations  can  predict 
displacement  profiles  which  are  qualitatively  correct.  It 
appears  that  it  will  not  be  difficult  to  predict  failure  by 
these  means.  The  AFTON-NONSAP  features  responsible  for 
inaccuracies  in  displacement  predictions  are  simplified  boun- 
dary conditions,  axial  symmetry  assumption,  and  premature 
termination  of  the  driving  pressure. 

3.  AFTON  results  can  be  simply  scaled  to  account 
for  variation  in  projectile  diameter. 
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SECTION  VI 

SUMMARY  AND  CONCLUSIONS 

Hydrodynamic  ram  effects  induced  by  high-velocity 
fragments  have  been  studied  with  experimental  and  analytical 
and  nun\erical  techniques.  The  projectiles  were  11.1  and  14.3  nm\ 
diameter  steel  spheres  traveling  in  excess  of  1.5  km/s. 

It  was  found  that: 

1.  High  velocity  projectiles  constitute  a severe 
threat  to  entrance  panels. 

2.  2024-T3  aluminum  alloy  is  far  more  resistant  to 
ram  damage  than  7075-T6  aluminum  alloy. 

3.  Both  ballistic  foam  and  styrofoam  are  extremely 
effective  in  defeating  entrance  panel  hydrodynamic  ram. 

4.  The  dynamics  of  entrance  panel  motion  and  failure 
are  as  follows:  The  impact  induced  shock  wave  produces  ex- 
tremely high  pressures  over  the  area  immediately  adjacent  to 
the  impact  site;  however,  the  duration  of  the  shock  pressure  is 
only  tens  of  microseconds.  During  the  time  that  significant  loads 
are  applied,  the  response  of  the  panels  is  largely  inertial.  Most 
of  the  impulse  is  transferred  to  the  panel  before  it  has  strained 
appreciably.  The  outward  motion  of  the  panel  then  proceeds  at  a 
relatively  slow  rate;  the  major  displacement  step  propagates  at 
only  'v-300  m/s.  The  hoop  stresses  due  to  the  deformation  induced 
by  the  shock  may  cause  very  early  tensile  failure  at  the  entrance 
site.  If  this  occurs,  cracks  propagate  out  as  the  panel  displaces. 
Thus,  the  critical  physical  processes  in  an  entrance  panel  event 
take  place  within  a few  tens  of  millimeters  of  the  entrance  site 
during  the  first  hundred  microseconds  after  impact. 

5.  The  local  nature  of  entrance  panel  phenomena 
permits  meaningful  experiments  to  be  carried  out  in  small  scale 
tanks . 
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6.  The  moire  fringe  device  developed  for  this  program 
can  be  used  to  measure  hydrodynamic  ram-induced  panel  motions 
with  sufficient  accuracy  and  resolution  to  check  numerical 
predictions . 

7.  The  fluid  drag  code  developed  at  NWC  provides  an 
adequate  description  of  fluid  pressure  in  the  tank  interior. 

8.  The  fluid  shock  code  embodies  assumptions  about 
wavefront  shape  and  velocity  decay  which  are  clearly  violated 
in  hydrodynamic  ram  events. 

9.  The  BR-lA  (HR)  code  is  not  able  to  adequately  predict 
hydrodynamic  ram-driven  panel  motions.  The  principal  reasons 
apparently  involve  the  basic  nature  of  the  piston  theory 
embodied  in  the  code. 

10.  When  driven  by  AFTON-generated  pressures,  the 
BR-lA  can  probably  provide  a sufficient  description  of  panel 
motions . 

11.  The  AFTON-NONSAP  calculations  of  panel  displacement 
are  qualitatively  correct.  Small  refinements  of  the  technique 
will  probably  yield  excellent  quantitative  predictions  of  panel 
shape  as  well  as  a failure  criterion. 

12.  2024-T3  aluminum  alloy  panels  1.5  mm  thick  fail 
catastrophically  when  struck  by  a 11.1  mm  diameter  sphere  above 
'■^1.5  km/s., 

13.  A layer  of  ballistic  foam  equal  to  the  projectile 
diameter  defeated  a 2.38  km/s  impact. 

14.  In  the  configuration  used  here,  no  side  panel  or 
exit  panel  failures  occurred. 

15.  Partial  results  indicate  foam  is  of  marginal 
utility  in  defeating  ram  effects  at  exit  panels. 
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APPENDIX  A 


This  Appendix  was  prepared  by  California  Research  and 
Technology  under  subcontracts  RI-76884  and  77713  for  the 
University  of  Dayton  Research  Institute,  Dayton,  Ohio. 
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FINITE  DIFFERENCE/FINITE  ELEMENT  ANALYSIS  OF  FUEL  TANK 
PENETRATION  BY  STEEL  FRAGMENTS 

Sheldon  Schuster 
Y.  Marvin  Ito 

California  Research  L Technology,  Inc.* 

Woodland  Hills,  CA  91367 

A. 1 INTRODUCTION  AND  APPROACH 
A. 1.1  Prior  Study 

In  a previous  study  by  CRT‘  , the  WAVE-L  two-dimensional 
(2-D)  Lagrangian  finite  difference  code  was  used  to  investigate 
the  dynamics  of  .50  caliber  bullets  penetrating  fuel  tanks, 
with  particular  attention  being  given  to  damage  to  the  exit 
panel,  i.e.  the  rear  surface  of  the  tank  which  is  penetrated 
after  the  projectile  passes  through  the  fuel.  The  impacts 
were  at  normal  incidence  with  no  yaw,  but  tumbling  was  simu- 
lated by  continuously  changing  the  shape  and  frontal  area  of 
the  bullet  to  match  an  assumed  time-dependent  drag  force. 

A . 1 . 2 Objectives  of  Present  Study 

In  the  current  computational  program,  a study  has  been 
performed  of  the  effects  of  steel  spheres  (representing  chunky 
fragments)  impacting  the  entrance  panel  of  aluminum  tanks 
(filled  with  water  simulating  fuel),  with  and  without  foam 
backing.  The  problems  which  were  numerically  analyzed  during 
this  study  are  shown  in  Table  A-1. 


* This  work  was  performed  by  CRT  under  Subcontract  RI-76884 
and  Hl-77713  with  University  of  Dayton  Research  Institute. 

‘ M.  Rost nblatt , G.  E.  Eggum,  and  L.  A.  DeAngelo,  "Numerical 
Analyses  of  Fuel  Tank  Penetration  Dynamics,  AFFDL-TR- 76- 31 , 
California  Research  ^ Technology,  Air  Force  Flight  Dynamics 
Lab April  1976 . 
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tart.f:  a-1.  CCWDITIONS  FOR  CRT  CALCOLAIIOSS  AliD 


The  objectives  of  tliese  calculations  were  (1)  to  imple- 
ment new  and  more  efficient  numerical  methods  for  analyzing 
fragment  and  projectile  impacts  into  fuel  tanks,  and  (2)  to 
validate  these  methods,  insofar  as  possible,  by  comparing 
numerical  results  with  dynamic  observations  of  the  penetration 
and  response  processes  occurring  in  tanks  impacted  by  steel 
spheres.  These  observations  were  obtained  by  UDRl  in  connec- 
tion with  the  experimental  program  described  in  the  main  body 
of  this  report. 

A.l..^  Numerical  Tecliniiiucs 


i 
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A two-phase  analysis  method  was  employed:  In  Phase  1, 
an  adaptation  of  the  Ab'TON  2-Ll  finite  difference  code*  was  used 
to  calculate  the  ponotvation  dunamioa  and  to  determine  the  time- 
resolved  force  loading  between  the  water  and  the  tank  wall.  In 
Phase  11,  the  tank  wall  veoponac  was  calculated  by  applying 
this  force  loading  at:  a boundary  condition  in  a NONSAP  2-D 
f inite -element  code*  solution.  This  approach  of  partially  de- 
coupling the  ponetra  on  and  response  assumes  that  tl\e  pene- 
tration dynamics  are  largely  unaffected  by  the  rigidity  of  the 
front  plate  (although  the  plate  mass  and  inertia  arc  correctly 
taken  into  consideration).  This  is  a reasonable  assumption 
during  tlio  early  stages  treated  by  the  Phase  I analysis. 

a.  Ain'ON  Finite  lUfference  Code 

TIk'  AFTON  code  was  selected  for  this  study  because  of 
its  ability  to  continuously  redefine  the  Lagrangian  mesh 
spacing.  This  capability  reduces  the  large  distortions  which 
occur  in  those  problems  when  the  mesh  moves  in  the  normal 
Lagrangian  manner  with  the  material.  In  the  Phase  1 calculations , 


* W.  J.  Niles,  J.  P.  Germroth,  and  S.  H.  Schuster,  "Numerical 
Studios  of  AFTON  2A  Code  Development  and  Applications,  Vol.  II," 
AFWL-TR- 70- 22  Vol.  II,  February  1971. 

* K-J  Bathe,  F..  L.  Wilson,  and  R.  H.  Iding,  "NONSAP  - A Struc- 
tural Analysis  Program  for  Static  and  Dynamic  Response  of 
Nonlinear  Systems",  Structural  F.ngineering  Laboratory,  Univer- 
sity of  California,  Rerkoley,  California,  DC  SF.SM  74-3,  Febru- 
ary 1974. 
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the  mesh  points  defining  both  the  front  surface  of  the  water 
and  the  cavity  behind  the  penetrator  moved  with  the  particle 
velocity  associated  with  the  material  (i.e.  Lagrangian) . 

However,  interior  mesh  points  (within  the  fluid)  wore  contin- 
uously rczoned  toward  positions  nominally  equidistant  from 
their  neighbors.  As  illustrated  by  the  plot  of  the  mesh  in 
Figure  A-1,  after  the  sphere  has  penetrated  about  one  diameter, 
the  severe  distortions  of  the  mesh  usually  associated  with  pure 
Lagrangian  calculations  are  minimized  by  the  AFTON  technique. 

At  the  some  time,  the  motion  of  the  boundaries,  whicli  would  be 
lost  if  the  mesh  were  fixed  in  space  and  time  as  in  pure 
Eulerian  formulations,  has  been  preserved. 

In  the  Phase  I calculations,  the  U.16-cm  (l/)6-in.) 
aluminum  entrance  panel  was  modeled  by  attaching  the  inertial 
mass  of  the  aluminum  to  the  mass  of  the  mesh  points  along  the 
front  surface  of  the  fluid.  Material  models  for  the  water  and 
foam  were  the  same  as  used  in  the  previous  study.'  The  steel 
sphere  was  treated  as  a rigid  body. 

The  Phase  1 calculations  were  stopped  before  the  first 
signal  reached  the  nearest  side  wall  of  the  tank  (at  r - 37.5  cm), 
and  when  pressures  everywlicre  along  the  entrance  panel  were 
less  than  10  bars  (-  150  psi).  The  pressure  vs  time  data  along 
the  front  panel  were  tlien  used  to  derive  the  force  loading  for 
the  subsequent  I’hase  11  NONSAP  structural  response  calculation. 

b.  NONSAP  Ivinite  F.lement  Code 

NONSAP  is  an  implicit,  3-D,  non-linear  finite  element 
code  which  was  selected  for  the  second  phase  of  this  study  in 
order  to  adequately  and  efficiently  model  the  large  deflections 
of  the  thin  entrance  panel.  (A  finite  difference  code  analysis 
of  this  thin  plate  response  would  be  difficult,  due  to  the  very 
small  time  step  imposed  by  the  explicit  stability  criterion.) 

The  finite  element  calculations  assumed  the  front  panel  was  a 
37.5-cm  radius  disc  of  2024-T3  aluminum,  rigidly  supported 
around  its  edge,  and  with  a circular  liole  in  the  center  tlic 
size  of  the  impacting  sphere.  The  zoning  and  material  model 
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a.  Mesh  Spacing 


Front  Surface  of  Water 
(Inertial  Mass  of  Aluminum 
Entrance  Panel  Added  to 
Surface  IIuoli) 


b.  Velocity  Vector  Field 


^Front  Surface 
of  Water 


scale ; 

1 ki  ./ecc 


Figure  A-1.  Mesh  Spacing  and  Velocity  Vector  Field  in  Water  at 
6.2  usee  After  Impact.  CRT  Case  A:  1.11-cm 
Steel  Sphere  Impact  at  1550  m/sec 
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used  are  illustrated  in  Figure  A-2.  No  failure  criterion  vnus 
specified  for  the  aluminum,  so  only  elastic  and  continuous 
plastic  deformation  are  allowed.  An  axisymmetric  model  of  the 
entrance  panel  was  assumed  (instead  of  the  trapezoidal • shaped 
panel  used  in  the  UDRI  experiments)  So  that  a 2-D  version  of 
NONSAP  could  be  employed.  This  allowed  better  representation 
of  the  impulsive  force  input  and  of  the  highly  nonlinear  re- 
sponse of  the  plate. 
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ENTR<\NCE  PANEL  GEOMETRY: 


2Q24-T3  ALUMINUM 
material  MODEL: 


rigid  support 
at  37.5  cm  radius 


F<r,t)  detannlned 
by  AFTON  analysis 
of  pressures  In  water 


CROSS-SECTIONAL  VIEW 

GRID  SIZE: 

No, Elements  No.  Nodes 

11.1  mm  56  283 

14.3  mm  55  278 

TIME  STEP: 

At  • 2 psec. 


Figure  A-2.  NONSAP  Finite  Element  Axisymmatric  Model  of  the  Entrance  Panel 
of  the  Fuel  Trnk 
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A. 2 RESULTS  OF  THE  CALCULATIONS  AND  EXPERIMENTAL  COMPARISONS 


A. 2.1  CRT  Case  A;  Impact  of  1.11-cm  dia  Steel  Sphere  at 

1S50  m/sec  into  Water  Tank  with  0.16-cm  (1/16- in.) 

Aluminum  Wall 

A. 2. 1.1  Calculations 

The  progress  of  the  1.11-cm  dia  steel  sphere  (5.5  gm) 
through  the  water  in  the  fuel  tank  at  6.2,  36.1,  and  78.2  nsec 
is  illustrated  in  the  series  of  velocity  vector  fields  in 
Figures  A-1,  A-3,  and  A-4,  respectively.  At  6.2  nsec  (Figure 
A-1),  the  sphere  has  i 'uetrated  almost  a full  diameter  and 
cavitation  has  begun.  As  the  sphere  continues  to  penetrate, 
more  and  more  mesh  points  are  automatically  activated  so  that 
by  136  nsec,  3500  points  were  included  in  the  problem.  Since 
the  motions  of  the  front  panel  were  of  primary  interest,  ail 
zones  at  depths  greater  than  10  cm  were  replaced  in  the  calcu- 
lation after  136  nsec  by  a reflecting  boundary  condition. 
(Reflection  from  this  boundary  would  not  bo  felt  at  the  front 
panel  during  the  time  of  interest.)  The  Phase  I calculation 
was  continued  to  178  ysoc.  By  this  time,  the  pressures 
along  the  front  panel  had  decayed  to  less  than  a bar  (*'  15  psi) 
everywhere  except  within  5 cm  of  the  axis,  whore  several  zones 
were  still  at  about  10  bars. 

Time  histories  of  the  instantaneous  peak  pressure  in  the 
water,  and  of  the  velocity  and  displacement  of  the  sphere,  are 
shown  in  Figure  A- 5.  After  impact,  the  calculated  peak  pressure 
in  the  water  (which  always  occurs  just  ahead  of  the  penetrating 
sphere)  decays  rapidly  as  the  initial  shock  wave  expands.  At 
times  greater  than  25  usee  the  peak  is  essentially  the  dynamic 
pressure  (1/2  pU^),  whore  U is  the  instantaneous  velocity  of 
the  sphere  and  p is  the  fluid  (water)  density.  The  expected 
dynamic  pressure  (determined  from  the  calculated  instantaneous 
velocity  of  the  sphere)  is  shown  by  the  smooth  curve  in  Figure 
A-5a.  As  is  seen,  the  peak  pressure  directly  calculated  by 
the  code  closely  follows  the  1/2  pU^  curve  after  25  psec. 
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Figure  A-4.  Velocity  Field  in  Wnter  at  78.2  iisec  for  CRT  Case  A: 
dia  Steel  Sphere  Impacting  at  1550  m/sec 
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VELeCITV'  (CM/MoECi  DISPLACE  (CM)  PRESSURE  (KSaRS) 

0.  100.  200  . 0.  10  . 20  . 0 . 20.  *10. 


The  oscillations  in  the  calculated  peak  pressure  corac  about 
because  the  pressure  periodically  decays  in  one  zone  and  in- 
creases in  the  next,  consistent  with  the  zone  engulfment  time 
of  the  pressure  wave  through  the  mesh. 

The  effective  drag  coefficient,  Cq,  of  the  penetrating 
sphere  can  be  estimated  by  fitting  the  calculated  velocity,  U, 
and  displacement,  D,  histories  in  Figure  A- 5 with  the  equation 

U - (A-1) 

where  ■ 1550  m/sec  (the  impact  velocity)  and  with  3 ■ .055. 
Cp  and  3 are  related  by 

Cm  - ■!  3 (A-2) 

^ 3 pw 

where  pj,  and  are  the  respective  densities  of  the  sphere  and 
fluid  (steel  and  water),  and  r^  is  the  sphere  radius. 

During  the  first  130  usee  after  impact,  the  above  fit 
leads  to  a value  of  Cjj  ■ 0.62.  Initially,  the  drag  coefficient 
is  much  higher  (Cq^  •“  1.2),  but  it  decreases  to  the  effective 
value  (Cj)  •’  0.6)  within  about  1 dia  of  penetration. 

Time  histories  of  the  pressure  5 cm  from  the  axis  of 
symmetry  and  at  2,5  and  10  cm  depth  are  plotted  in  Figure  A- 6. 
The  duration  of  the  first  pulse  broadens  with  depth  as  the  time 
difference  between  the  arrival  of  the  initial  wave  and  the 
relief  waves  originating  from  the  front  surface  increases. 
Figure  A- 7 shows  peak  pressure  vs  radius  points  within  the  tank 
and  along  the  front  surface;  these  decrease  with  distance  from 
the  axis  for  each  depth. 

Pressure- time  histories  along  the  entrance  panel  out  to 
178  usee  from  the  Phase  I AFTON  finite  difference  calculation 
were  used  in  Phase  II  as  the  forcing  functions  in  the  NONSAP 
finite  element  code.  These  forces  were  applied  to  the  inside 
surface  of  a 37.5-cm  dia  disc  of  0.16-cm  (1/16-in.)  thick 
aluminum,  as  indicated  in  Figure  A-2.  The  calculated  displace- 
ment profiles  of  the  aluminum  plate  at  200  usee  intervals  are 
plotted  in  Figures  A-8  and  A-9.  The  peak  displacements  ranged 
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Tine  Histories  of  the  Pressure  In  the  Tank.  CRT  Case  A:  1.11-cm 
Dla  Steel  Sphere  Impacting  at  Vq  • 1550  m/sec  at  Points  5 cm  off 
the  Axis  of  Symmetry  and  at  2,5  and  10  cm  Depth  from  the  Entrance 
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7.  Peak  Pressure  Experienced  along  the  Front 
Panel  and  at  Selected  Points  in  the  Tank 
as  a Function  of  Distance  from  the  Axis 
of  the  Spherical  Fragment.  CRT  Case  A: 
1.11-cm  dia  Steel  Sphere  Impacting  at 
1550  m/sec 
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Figure  A— S-  Calculated  Dlsplacenent  Profiles  along  the  0.16  cm  Aluminum  Entrance  Panel 
at  r.equence  of  Tlir.es  Between  200  and  1000  jsec  After  Impact.  CRT  Case  A: 
1.11-cm  Dia  Steel  Sphere  Impa-'Clng  at  Vq  = 1550  m/sec 


from  about  2.2-cm  near  the  entry  hole  to  about  1-cm  at  radii 
between  10  and  30  cm. 

The  calculated  tensile  hoop  stresses  which  developed  in 
the  impacted  aluminum  plate  as  it  deforms  are  sliown  in  Figure 
A- 10.  Close  to  the  axis,  these  stresses  rcacii  a maximum  quite 
early  (t  < 200  psec)  and  then  decay  as  more  of  the  disc  is 
engulfed  by  the  shock  wave.  Only  a small  region  (out  to  about 
4-ciii  in  radius)  of  the  disc  attains  stresses  above  the  nominal 
70  kpsi  tensile  limit  of  2024-T3  aluminum.  It  is  tliereforc, 
reasonable  to  expect  that  any  radial  cracks  which  might  form 
around  tl\c  hole  will  not  propagate.  It  is  pointed  out  that 
the  peak  stresses  and  strains  which  would  lead  to  possible 
failure  develop  early  in  the  analysis,  prior  to  the  end  of  the 
Phase  I calculation. 

A.  2. 1.2  Experimental  Comparisons 

The  conditions  for  UDRI  Test  ET5B  correspond  approximately 
with  CUT  Case  A.  Experimental  observations  of  the  dynamic 
distortion  of  the  entrance  panel  were  made  in  this  test,  using 
the  experimental  arrangement  in  the  sketch: 


Entrance  panel  observed 
by  Moire  interference 
method  in 

-vertical  plane- — ^ 

-horizontal  plane>v.,^^ 


Sphere  impacting 
on  horizontal  O’ 
trajectory  — 


Figures  A- 11  and  A- 12  show  sequences  of  entrance  panel 
profiles  in  horizontal  and  vertical  planes,  as  observed  by  the 
Moire  interference  technique.  The  nominal  inter- frame  time  in 
these  observations  was  160  ysec.  The  impact  time  is  thus  known 
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Figure  A-12.  Experimental  Displacement  Profiles  Measured  Along  a Vertical  Line  through 
the  Impact  Point.  UDRI  Test  FT5B:  1.11— cm  Dia  Steel  Sphere  Impacting 
at  1550  m/sec 


to  within  ±80  nsec,  so  there  is  also  a ±80  nsec  time  uncertainty 
associated  with  each  profile.  The  displacements  in  the  vertical 
plane  are  consistently  larger  than  those  in  the  horizontal 
plane,  a result  of  the  plate  geometry  and/or  some  uncertainty 
in  the  displacement  observations. 

Taking  these  uncertainties  into  account,  Figure  A-13 
shows  that  the  calculated  profile  at  600  nsec  is  quite  similar 
to  the  observed  profiles  at  570  ±80  nsec.  However,  as  seen  in 
Figures  A-11  and  A-12,  the  observed  profile  continues  to  move 
out  after  600  nsec,  while  the  calculated  motion  (Figures  A-8 
and  A“9)  begins  to  slowly  rebound.  Assumptions  in  the  calcu- 
lation, and  differences  between  the  experimental  and  calculated 
conditions,  are  responsible  for  the  divergence  after  600  nsec: 

0 The  forcing  function  used  in  the  Phase  II  response 
calculation  was  set  to  zero  after  178  nsec,  the 
end  of  the  Phase  I analysis.  At  later  times,  there- 
fore, the  calculated  tensile  stresses  in  the  aluminum 
entrance  panel  were  able  to  arrest  the  outward  motion 
of  the  panel  and  turn  it  around  more  quickly  than  in 
the  experimental  case.  While  the  calculated  pressures 
at  the  surface  of  the  water  had  dropped  to  the  1 bar 
(-  15  psi)  level  or  less  by  178  nsec,  the  water  was 
still  moving  outward.  The  entrance  panel  in  the 
tests  absorbed  this  momentum,  thereby  reaching  larger 
displacements  before  rebounding. 

0 The  calculations  treated  the  panel  as  a flat  circular 
disc,  rigidly  clamped  at  the  outer  edge.  The  experi- 
ment used  a trapezoidal  entrance  panel  with  edges 
attached  by  bolts  to  the  side  walls  of  the  tank. 

This  attachment  presumably  strongly  restrains  edge 
motion  mormal  to  the  front  panel  surface,  but  may 
only  weakly  restrain  rotations  of  the  the  edge  or 
motions  in  the  plane  of  the  entrance  panel  (i.e., 
"radial”  motions).  This  relatively  weak  boundary 
constraint  would  allow  larger  displacements  in  the 
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test  case.  For  example,  based  on  large  deflection 
plate  theory** , the  maximum  displacement  would  be 
expected  to  increase  by  at  least  a factor  of  two 
if  the  boundary  is  allowed  to  move  freely  in  the 
radial  direction.  The  actual  boundary  conditions  on 
the  test  panel  were  somewhere  between  the  extreme 
of  the  immovable  boundary  assumed  in  the  calcu- 
lation and  a boundary  allowing  free  rotation  and 
radial  motion.  The  effect  of  the  geometry  and 
boundary  restraint  is  further  seen  in  a static 
solution  of  the  plate  deflection  in  which  the  0.5 
psi  hydrostatic  load  produced  by  the  water  head  at 
the  level  of  the  impact  point  was  uniformly  applied 
to  the  radially-restrained  disc.  This  solution 
produced  a O.S-cm  maximum  deflection  in  the  plate. 

By  contrast,  the  maximum  static  deflection  of  the 
test  plate,  as  measured  by  UDRI  with  Moire  finge 
shifts,  was  about  1 cm,  indicating  that  the  trape- 
zoidal shape  and  the  relatively  weaker  radial 
restraint  in  the  test  allows  significantly  larger 
deflections . 

The  entrance  panel  in  UDRI  Test  FT5B  showed  permanent 
distortion,  but  radial  cracks  did  not  propagate  outward  from 
the  entrance  hole.  This  appears  consistent  with  the  hoop 
stresses  calculated  in  CRT  Case  A.  The  plots  in  Figure  A-10 
show  that  calculated  tensile  hoop  stresses  exceeded  the  nominal 
ultimate  strength  of  2024-T3  only  near  the  entrance  hole. 


S.  Timoshenko  and  S.  Woinowsky-Krieger , "Theory  of  Plates 
and  Shells,"  2nd  Edition,  p.  404-415,  McGraw-Hill,  New 
York,  N.  Y.,  1959. 


A. 2. 2 CRT  Case  B;  Impact  of  1.11-cm  dia  Steel  Sphere  at 

1550  m/sec  into  Tank  with  0.16~cm  Aluminum  VJall 

Backed  by  3. 8 -cm  of  Foam 

This  case  was  identical  to  CRT  Case  A except  for  the 
inclusion  of  a 3.8-cin  layer  of  AVCO  ballistic  foam  placed 
directly  behind  the  entrance  panel.  This  foam  has  a crushing 
strength  of  about  25  psi. 

The  AFTON  finite  difference  calculation  was  begun  with 
the  steel  sphere  at  the  foam-water  interface  (i.e.,  it  was 
assumed  that  no  energy  or  momentum  was  lost  in  passing  through 
the  foam) . The  inertial  mass  of  the  aluminum  entrance  panel 
was  added  to  the  outermost  foam  zone.  At  a time  of  77  ysec. 
Figure  A-14,  the  flow  field  in  the  water  is  quite  similar  to 
that  in  Case  A (Figure  A-4).  The  entrance  panel,  being  separ- 
ated from  the  water  by  the  foam,  has  not  yet  felt  any  signi- 
ficant pressure  and  has  not  moved. 

By  174  ysec,  the  maximum  pressure  on  the  aluminum  is 
30  psi  and  only  the  foam  out  to  8-cm  radius  has  experienced 
pressures  sufficient  to  cause  any  crushing  (i.e.  P > 25  psi). 

In  Figure  A-14,  the  velocities  of  the  water  near  the  entrance 
are  such  that  the  water  appears  to  be  cutting  through  the 
first  few  centimeters  of  foam.  This  is  consistent  with  the 
observation  in  UDRI  Test  FT26B  that  about  1-cm  of  foam  has  been 
chipped  away  from  around  the  impact  hole. 

The  much  lower  pressures  generated  along  the  foam- 
backed  plate  surface  in  CRT  Case  B,  as  compared  with  the 
pressure  on  the  bare  plate  in  CRT  Case  A,  will  clearly 
produce  lower  displacements  and  stresses  in  the  aluminum. 

Since  the  calculation  of  the  no-foam  case  (Case  A)  showed  only 
a minimal  chance  of  crack  formation  around  the  entry  hole, 
and  since  examination  of  the  post-test  entrance  panel  from 
the  corresponding  UDRI  Tests  (FT26  and  FT26A)  showed  no  mea- 
surable permanent  distortion  (other  than  the  entry  hole) , it 
was  decided  that  a phase  II  NONSAP  calculation  of  Case  B would 
net  be  useful . 
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Figure  A-14.  Velocity  Field  at  77  psec. 
CRT  Case  B:  1.11-cm  Dia  Sphere  Impacting 
* * * at  1550  m/sec.  The  Water  is 
«####.  . Separated  from  the  Entrance 
Panel  by  3.8-cm  of  AVCO 
tit,,,,..  Ballistic  Foam 
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A. 2. 3.1  Calculations 

The  conditions  for  this  case  are  identical  to  CRT  Case  A 
except  that  the  impacting  sphere  diameter  is  1.43-cm  instead  of 
1.11-cm.  If  the  thiokneae  of  the  entranoe  panel  is  diavegavdedt 
the  penetration  dynamics  in  Case  C are  seen  to  be  identical  to 
those  in  Case  A at  linearly  scaled  times  and  dimensions.  The 
scale  factor  is  the  ratio  of  the  diameter  of  the  spheres,  i.e., 
1.43/1.11  = 1.29.  Thus,  for  example,  the  pressure  and  particle 
velocity  at  2-cm  depth  at  10  ysec,  as  calculated  for  the  CRT 
Case  A conditions,  would  be  the  same  as  those  found  in  the 
Case  C impact  at  2.6-cm  depth  at  12.9  psec. 

These  considerations  suggest  that  it  may  be  possible 
to  use  scaled  results  of  Phase  I penetration  analyses  obtained 
for  one  size  fragment  or  sphere  to  drive  Phase  11  tank  wall 
response  analyses  for  other  size  fragments  or  spheres.  This 
approach  could  be  attractive  for  future  design  studies,  so  it 
was  evaluated  here  by  making  two  analyses  of  Case  C.  In  the 
first,  which  will  be  referred  to  as  CRT  Case  C-1,  appropriately 
scaled  Phase  I penetration  and  entrance  panel  loading  results 
from  the  CRT  Case  A calculation  were  used  to  drive  the  Phase  II 
NONSAP  analysis  of  the  entrance  panel  response  to  the  larger 
sphere  impact.  In  the  second  analysis,  referred  to  as  CRT 
Case  C-2,  an  AFTON  Phase  I analysis  of  the  penetration  of  the 

1.43- cm  dia  sphere  was  obtained  (using  the  correct  inertial 
mass  for  the  0.16-cm  thick  aluminum  entrance  panel).  The 
validity  of  the  scaled-loading  approach  used  in  Case  C-1  can 
be  evaluated  by  comparing  the  results  of  C-1  and  C-2. 

a.  Case  C-1:  Analysis  Using  Scaled  Loading 

The  entrance  panel  thickness  is  the  same  in  CRT  Cases 
A and  C,  while  the  sphere  diameter  increases  from  1.11  to 

1.43- cm.  The  relatively  thinner  plate  (with  respect  to  the 
sphere  diameter)  in  Case  C thus  provides  less  inertial  con- 
finement to  the  water  during  the  penetration.  As  a 
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consequence,  the  scaled- loading  approach  results  in  some  over- 
loading of  the  plate  in  Case  C-1. 

Displacement  profiles  obtained  at  200  psec  intervals  from 
the  NONSAP  analysis  of  the  entrance  plate  response  are  shown 
in  Figure  A- 15.  These  have  the  same  general  shape  as  those  for 
the  smaller  sphere  in  Case  A,  (Figure  A-8)  but  displacements 
are  consistently  50t  greater  in  magnitude;  3.5  vs  2.25-cm  on 
axis,  1.5  vs  1.0-cm  at  15-cm  off  axis,  etc.  Similarly  stress 
profiles  (Figure  A-16)  shows  consistently  higher  stresses 
relative  to  those  at  200  psec,  for  example,  about  2.5  times  as 
much  area  of  the  front  panel  (5. 75  vs  3.75-cm  in  radius)  was 
stressed  by  hoop  tension  greater  than  the  nominal  70  kpsi 
tensile  limit  of  the  2024-T3  aluminum  than  in  the  previous  case 
(Figure  A-12).  Also,  at  400  Psec  the  hoop  stress  level  i. 
still  just  below  the  tensile  strength  as  far  out  as  10-cm  radius, 
which  shows  a greater  tendency  for  crack  propagation  than  the 
previous  case.  Perhaps  the  most  definitive  change  occurs  in 
the  hoop  strain.  As  shown  in  Figure  A- 17,  the  strains  produced 
at  200  psec  by  the  scaled  forces  in  Case  C-1  are  considerably 
larger  than  those  in  Case  A.  In  the  first  calculation  the 
strain  drops  rapidly  from  a peak  of  91  to  under  2%  at  1.2-cm 
and  is  less  than  0.5%  by  3.5-cm  while  the  larger  sphere  pro- 
duced a peak  strain  over  13%  and  the  2%  and  0.5%  levels  occur 
at  almost  3 and  7-cm  respectively.  Thus  at  least  six  times  as 
much  material  underwent  hoop  strains  of  2%  or  more  in  the 
scaled  calculation,  a difference  great  enough  to  provide  a 
potential  failure  criterion. 

b.  Case  C-2:  Analysis  Using  Loading  Calculated  for 
1.43-cm  dia  Steel  Sphere  Impact 

A Phase  I AFTON  code  analysis  was  made  of  a 1.43-cm  dia 
steel  sphere  impact  at  1500  m/sec.  (This  slightly  lower  impact 
velocity  corresponds  more  closely  with  conditions  of  UDRI  Test 
FTA9,  which  will  be  used  for  comparison  later.)  Comparisons 
of  peak  pressure.  Figure  A-18a,  and  the  specific  impulse  up 
to  130  psec,  Figure  A-18b,  for  Cases  A,  C-1,  and  C-2  indicate 
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Distance  Irom  Point  o£  Impact  (cm) 

Figure  A-'IB.  Peak  Pressure  and  Specific  Impulse  to  130  ysec  along  the 
Entrance  Panel  Calculated  by  the  AFTON  Code  for  CRT  Cases; 
A (l.ll>cm  Dia  Sphere),  C'*!  (Time  and  Distance  in  Case  A 
Scaled  by  1.43/1.11),  and  C-2  (i.A3'~cm  Dia  Sphere) 
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that  simple  scaling  of  the  forces  from  the  Al-TON  calculation 
in  Case  A provided  a reasonable  approximation  to  the  actual 
forces  calculated  in  Case  C-2,  Initially,  the  cavity  behind 
the  penotrator  appears  to  bo  the  primary  source  of  relief  waves. 
The  scaled  peak  pressure  lies  along  that  calculated  for  the 
1.43- cm  diameter  sphere,  while  the  scaled  integrated  specific 
impulse  is  about  201  higher.  At  distances  large  relative  to 
the  siae  of  the  ponetrator  (about  8-cm)  the  relief  from  the 
entrance  wall  dominates  the  attenuation  and  the  peak  pressure 
is  almost  the  same  for  both  spheres.  Similarly  the  difference 
between  tlio  scaled  and  actual  impulses  grows  with  range  so 
that  at  10- cm  the  scaled  value  is  almost  a factor  of  two  too 
high.  Since  the  stresses  and  strains  loading  to  fracture  occur 
early  in  time  and  close  to  the  entrance  hole,  results  of  Case 
C-1  are  probably  adequate  to  describe  the  response  to  the 
1.43-cm  sphere.  The  calculated  displacements,  liowever,  are  a 
result  of  long  time  and  largo  range  pressures,  and  are  probably 
higher  than  would  bo  expected  from  tlie  C-2  forces. 

A.  2.3. 2 lixper imentai  Comparison 

Observations  of  entrance  panel  dynamic  deformation  for 
conditions  corresponding  to  CRT  Case  C were  not  obtained  in 
the  UDRI  tost  series.  However,  measurements  were  made  in  UDRl 
Test  FTA9  of  the  displacement  of  the  sphere  vs  time.  Figure 
A- 19  compares  the  measured  displacements  (which  wore  obtained 
to  about  55  usee)  , with  the  displacements  calculated  in  CRT 
Case  C-2.  Excellent  agroomont  is  seen;  the  slight  offset  in 
time  is  probably  due  to  the  difficulty  in  establishing  the 
exact  zero  time  by  forward  extrapolation  of  the  expurimental 
data. 

Shock  front  arrival  times  at  increasing  radii  along  the 
entrance  panel  were  also  measured  in  UDRl  Test  FTA9.  These  are 
compared  in  Figure  A- 20  with  the  times  at  which  peak  pressures 
occured,  as  calculated  in  CRT  Case  C-2.  The  length  of  the 
vertical  bar  associated  with  each  calculated  peak  corresponds 
to  the  width  of  the  corresponding  zone  and  is  a measure  of  the 
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Position  (cm) 


Distance  from  Point  of  Intact  (cm) 


CRT  Case  C-2:  1.43-cm 
dla  Steel  Sphere  Impact 
at  1500  m/sec V 


t’DRI  Test  FT9A:  1.43-cm 
dla  Steel  Sphere  Impact  at 
1460  m/sec 


Time  (visec) 

Figure  A-20.  Arrival  Time  of  Experimental  Shock  Front  and  Calculated 
Peak  Pressure  along  Entrance  Panel 


! 


uncertainty  in  position.  Again  the  agreement  is  seen  to  be 
excellent.  Since  shock  velocity  and  hence  displacement  of  the 
shock  front  along  the  front  panel  is  a function  of  pressure, 
this  agreement  indicates  that  the  calculated  peak  pressures 
are  correct. 


A. 3 SUMMARY  AND  CONCLUSIONS 


1.  A two-pha  e numerical  technique  for  predicting  pene- 
tration and  response  of  fuel  tanks  impacted  by  high  velocity 
fragments  has  been  demonstrated  and  correlated  with  available 
experimental  observations.  Phase  I of  the  method  uses  an 
adaptation  of  the  AFTON  2-D  finite  difference  code  to  analyze 
penetration  dynamics  and  to  determine  the  time-resolved  pres- 
sure applied  by  the  fuel  along  the  entrance  panel.  Phase  II 
uses  this  information  to  drive  a NONSAP  2-D  finite  element 
code  analysis  of  the  entrance  panel  dynamic  response  (defor- 
mation) . 

2.  Using  this  method,  calculations  have  been  made  of 
steel  sphere  impacts  at  ^/1500  m/sec  into  tanks  containing  water. 
Results  of  the  calculations  vrere  compared  with  UDRI  experi- 
mental observations  of  projectile  displacement  history  within 
the  water,  shock  arrival  time  in  the  water  along  the  water- 
entrance  panel  interface,  and  entrance  panel  displacement 
profiles.  Excellent  agreement  was  obtained  between  the  calcu- 
lated and  experimental  projectile  motions,  as  well  as  the  shock 
arrival  times.  Comparisons  of  entrance  panel  deflection 
profiles  are  generally  good  at  early  times,  but  the  measured 
profiles  eventually  show  larger  displacements  than  the  calculated 
profiles.  This  reflects  differences  in  boundary  conditions 
between  the  tests  and  calculations,  as  well  as  the  fact  that 
the  penetration  analysis  (Phase  I)  was  not  continued  far 
enough.  (This  is  not  a limitation  of  the  technique;  AFTON 
solutions  of  penetration  could  be  extended  in  future  analyses 

if  desired.) 

3.  Catastrophic  failure  of  entrance  panels  occurs  by 
propagation  of  radial  cracks  emanating  from  the  impact  hole. 

The  calculations  show  that  peak  tensile  hoop  stresses  and 
strains,  which  would  lead  to  such  rupture,  develop  in  the 
region  around  the  impact  hole  well  aitzK  the  impact,  but  well 
bz^oKz  most  of  the  panel  responds  to  the  impact.  In  the 
problems  considered  in  the  current  program  (1.11-cm  and  1,43-cm 
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dia  sphere  impacts  at  ''1500  m/sec). 


these  peak  stresses 


devel- 


oped within  200  psec  of  the  impact. 


4.  Use  of  foam  backing  on  the  entrance  panel  greatly 
reduces  the  transient  pressures  applied  by  the  water  to  tt.e 
entrance  panel. 
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RECOMMENDATIONS 


A.  4 

This  numerical  approach  should  be  used  to  analyze  inter- 
actions and  mechanisms  which  are  difficult  to  identify  or 
observe  experimentally,  and  which  may  be  important  in  devel- 
oping tank  designs  which  are  more  resistant  to  impact  damage. 

An  important  example  would  be  to  use  numerical  solutions 
to  systematically  examine  the  effects  of  foam-backing  para- 
meters on  entrance  panel  response.  Foam  backing  is  clearly 
very  effective  in  reducing  damage,  but  its  use  reduces  the 
available  tank  volume.  For  this  reason,  it  is  highly  desirable 
to  optimize  the  foam  system  by  choosing  the  best  combination  of 
foam  thickness,  foam  crushing  strength  and  foam  density. 

Systems  in  which  foam  or  honeycomb  are  sandwiched  between 
layers  of  aluminum  (or  other  materials)  also  can  be  evaluated 
and  optimized  in  this  manner. 


APPENDIX  B 


This  Appendix  provides  a listing  of  the  changes  made  to 
I the  PDC  and  the  listing  for  the  shock  propagation  code.  None 

of  the  changes  in  the  FDC  affected  the  way  the  code  calculates 
pressure.  They  were  made  to  adapt  the  code  to  the  UDRI  tank 
boundary  conditions,  and  to  provide  output  in  the  format 
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Figure  B-1.  Listing  of  modified  portions  of  subroutine 
IMAGES . 


! 

i 

j 

'2 

•i 


IBIS  TASS  IS  BEST  QUALITY  PRACTlCABSj| 
IBQM(»^rXE,WISHEiD10DD,Q 


1 P(!nr.»A“  HtJC*  <TN“'IT, OUTPUT, T 40.  ««I»Jo()T,T4P'"(i,0()'jpuT,x(ip- 

..aiMtaSlON 

'’lUENSIOU  KOT(«>t,ifT(‘.) 

CQMMQNIPXQI,ASXiW.LJ.aOi-*ltal.*^UUYCt5JiJ  .TilTOLl  ,X1  (3:H  .Yi  UJ.l  .2H 

S lU')  .Ur.,OT,TM4Y,Nto,yMtT3;.M  .BOEl('-o.JC:).»'’I‘?ST(5:) 

CnxHON  /HQ£L/  yjQn.NJQO 

oo-MOM/iniKira/ipif,) 

’.'LAfllT'iiaOD.itJQa  . - ... 

00  1 YYil.N  100 

an .XJo.!l  •?  XY 

PF40  (F.PaCUUP 

»!F4P(g.30gllTHty.nT.TPLnT , 

3]ei  •'Ot  MAT(2Fi;.r,,iij 

ilQ21*J.Aii?  . 

15  PfACK?,  too  II  (Xwei.j)  , J»l,  3> 

2- -ilQUTiNUi .. 

tool  F0^H4T(3F1 J.5> 

gfloi  FQitMOTtlTI 

1000  FryBM4T(Il) 

.AJl ...C.ALL...luiGt.(N3.t._...  

04LL  T04J(  4,0,0*,y;.,N  ,r.,n  US,  JYi,.l'44Yt 

OALL  JXl.  JMAX) 

1 OONTTUUS 

P^INT*.  -riATA  TO  .PiLNOM" 

25  ITM4y»  TH4y»10',  ♦ .JOOOl 

,0.0  .n  rriMi  ..•.  iatMAx  . . 

IFOUUT  * 0 

_P.0.15  .I“OINr«t.HP  ... 

TF  ( IFnsTdPflNTI  .UT,  rTIS-  ) tCOUNT  i TCOIINT  ♦ I 

ju 15 oomilUi 

TIMt  ■ ITTU-r»10 

. .,.«RI.TtlT,l?,l  Tl‘'E.,,l.conNT 

HRIT-(h,lUl  TI“i,ICOIlNT,ITIMi 

i;.  .'■p''H4T(.vt';,AiI''J. 

35  U F0FM4T<1X,  f1q,(,,I«;,.30X,If,1 

X > ^ 

L»0 

IF  .(10 PUNT.  vO. A t GOTO  Xl 

IFLAOsO 

All !10.A5_MaiJ!iUi! _ . . 

YOTtlUOi-yi  s c 

XTdNOlXl  » J 

35  OONTINUF 

O.Q_:»A_IMX!it  - . . 

*5  TF  ( TFlAO,  ;n,  i ,4Nn,  IPOINT,  fo.NO  ,ANn,  IFtojithfOINI  » ,r,  ^.itihh 

__1 opto  21 _ 

TF  (ifirftitpoiuti  ,r.:.  itthf)  r.n  to  to 

IF  tL.£Q.C  .OR.  L.LO.f.  .OR.  L.CO.ig  .OF.  L.Fp.lj  .QR.  L.FQ.gQ 

1 .OR,  L.M.TO  .0».  L.E0.35  ,0f,  L.tO.U-*  ,0»,  L.tO.ltl 

2 L.J(_!_3 

IFLAG  . t 

<t~*  1 ' 

XPTKl  . TPOIMT 

55  KTIXI  I ITIMF 

IS...  .L -X-a.Ift.  T0_Z1 ..  . 

TF  (IPOrn.-O.NPI  GOTO  21 
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program  HRAM  and  MIRROR. 
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IF  ( KQUNT  .EO.  5 ) KNT  ■ 2 

IF  i KOUNT  .EO.  8 1 KNT  t S'  . - 

« IF  ( KMT  .CT.  0 t CALL  PICTURE 

CALL  PICTURE  - ' 

IF  ( NS  .LT.  l.L  > GO  TO  S 

ST«P 

END 


FUNCTION  RHOl  (A«NI 
COHHJN  /If  AINCtBINC 
EXTERNAL  HSF 
REAL  NSF«N 
A1  « AINC 

A1  --•-ftlNft 


0 • 1,  ♦ lE-3 
10  C • NSFtaiNI 

0 ■ A • C 

IF  < 0 .LT.  1£-12  .A.  0 

* PRINT  TiOO.  A,  B.  C.  D 

-»« FOR(«rH“-A»' 

IF  I 0 .LT.  0.  I 

IF  < 0 .GE.  Al-  » 

IF  ( Al  ,GE.  lE-13  » 

PRINT  1000  . 0.  Al  - 

GO  TO  80 

M - IF-R  ■Rl-.Gt:.— 1-E-ig  > -GO  TO  31 

PRINT  1001,  0,  31 

- GO  TO  89  - 

30  9 « B - 91 

ei  • 31  /-  2r  

GO  TO  NO 


6T.  -lE-12  1 


GO  TO  60 


GO  TO  20 
GO  TO  SO 
GO  TO  60 


>0  3 » 9 ♦ 31 

GO-  TO  10  - - - - 

60  CONTINUE 

RHOl  • -i  — 

RETURN 

• FRON  POINT  1 *•  > 

1001  PORNAT  <-  ERROR  RETURN  NiSSAGE*  •-"EIE.S'’ 

• FROH  POINT  2.-  I 

- - EOT 


-At  •••ei2.-5*'-«ETURN€»- 


, B1  p-ElE.B-  RETURHEO 


Figure  B-4.  Listing  of  Shock  Propagation  Code. 
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fflis  Mm  IS  BSSI  asialK  mO«»JM 

jjgOjl  COP^  j'UBlNlSHEi)  TO  DDO 


SUnKOUTtNE  S<STUP 

^nOHHON^^l/  KTAMS,0«H,TH,BTM,N.MS,Pt,9EU.U,ttH.mJOO>,OUT,OOHH 

COiNON  AJ/  AXNCtBZNC 
REAL  HS.N 

6*  T*-S+ AM5-IN-eotyHH9-i»>-e»,  - 3Ss-AHO--A»s 

DATA  BETA  ,AINC  ,eiNC  .LI  / 

"*^®***  »!••  / 

DATA  H ,HS  / 

• r.  .1.  / 

• R«I»  HAS  LI  VALUES  IH  THE  IHTERWAL  < 1,  0 I 

HAHTHHtA-TO  BE  t. 

LtPt  ■ U ♦ 1 

Rill  ■ It  _ 

no  10  I ■ 2,  Ll 

RIII  • FLOATtLiPl  - 11  / FLOATALU 
10  CONTINUE 


PUNCTIO'4  UlUt 
ui  ■ It  • It/A 
RETURN 
ENP 


REAL  fUNOTlOS  HSF  < Bt  X t 

• - • HSF  » I n»*<x»i.i  '-  0 I -/  < x«(e-it>  1 

HSF  B sa<T(HSF> 

RETURN 

END 


FUNCTIPH  DHHOHRH) 

COHHON  /!/  eETAHSiORH«THiBTHtNtHStPliaETAiLl>tlRliR($IIO)  tDUTtDDRH 

• ,KNT 
REAL  HStN 
RHN  B RH  **  N 

— ^f-*-N-B-HS-  *-  -TPH- 1 ( ♦**« — 

HTM  • N • RHN  • RH  - (N»l.>  • RHN  ♦ It 

ORHOl  B TOP  / BTH 

RETURN 

ENO  - 

FUNCTION  OUTHTRHI 

COHHON  /!/-  KTA.HStORHtTHtBTHtNtHStPltOETAtlltLlPl.Rt300»,OUT,OORH 

• tKNT 

- REAL-NS-fH ■ 

OORH  B oorhokrhi 

-- OUTN-b-  t 'JORH-«-tt)RWeRM»  FRH  t -»-BRW 

nUTH  B OUTH  * (META/RHI 

• RETURN - ■ ■ • 

END 


Figure  B-4  (Cont'd) . Listing  of  Shock  Propagation  Code, 
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IBIS  PAGQS  IS  BUST  QUALITY  mCIlCABj^ 
I5j(0|BCI0^PYP;^RiMSB^»TO^;^Q  ' 


PROGRAM  VKPRSSS  <INPUT«OUTPUTtPLOTI 

*XT€RNAL  MSP  

DIMEliSlON  REQtSOOt  ,T2e«3aOI.T3PI300).Tl>Bt3SOt 

COiNON  /!/  B:rAMS«aRH.lH«ilTH.N.MS.Pl«8£rA.LltllPliR(30Q>«DUTtOORH 

• ,<NT 

— COHMOM-/2y  «><300+-»«^S(W-fW2««lrt-*>2«***vT21<3tM-»T2f+a««1 

• ,T3C300>  «T<>(300) 

DIMEMSIO')  XMSISOt  - . . 

REAL  N.NStMSP 

DATA  XMS/1.09S>t«0&«l.a34it>029<1.02ttt«01A,l>Oi3<l>OOS/ — 

CAIL  SETUP 


KOUNT  - KOUNT  ♦ 1 

MS  ■ XM5(KOUNT»  - - 

IF  t MS  .£0.  0.  I STOP“MS  • 0." 

DEFINE  THE  EXPRESSION -F«R  Tli-#-Tl--«-Pli 

RH  « RHOKHSiNI 

U ■ UURHt — 

ORH  « nRHOKRM.MSiNI 

TH  ■ SETA  * MS  • ORH 

TH  ■ TH  / <U»RHli 

QP2  I 3.»RH  - 1. 

SOUVE  FOR  THE  TIRHSI  P2,  P3t  AND  PM  OF  EQ.  70 

t>S  ■ U « -U " 

UTH  « U • TH 

USTH  ■ -US  ♦ Trt 

URLTA  « U • BETA 

USTH3  • UTH  ♦ UTH- - ■ 

UTHB'ETA  • UTH  ' BETA 

-OOHHle  «-«-8a-T^  * MS  «-BHH-F  F t RH  » RH  t 

P2l  « US  - U ♦ UBETA 

- P22  • PH  * CONNIE- - 

P2  » P22  ♦ UTH  • USTH 

• 8ETAHS-*-«£TA-«-MS 

OUT  « nUTHlRHJ 

P3  ».  UTH  t USTH5  - imUSTH  --UTHBETA— -B&TAMC  

PM  » U5THS 

find  THE  TERMS  T2(I»»  TJ«I>,  ANO-TAUI 

T2A  * PH  / Q<»2 

T3A  ■ -TZA  /-0P2 

TMA  » tJA  ! QP2 

-BO  -AB  I - t 

REQ(H  « R(It  ••  QP2 

T28«I»  • i.—  -REOdi — 

X < Rdl 

• CHPIS  • OP2  ♦ -AEOCTXT - 

CMFI3  « QP2  • AlOGtXi  - 1. 

T.3E44F-  *-4-.  PC04  AF  -P-CMR-IS 

TMBII)  ■ 2,  - REad»  * ( CHRIS  » CHRIS  » 1.  I 

T2»I)  » T2A  • P2  • T2Bd» — 

T21(I»  * T2A  * P21  • T2Bdl 

T22<I»  TZA  • P22  « T2»d»  - - 

T3d»  • T3A  * P3  • T3«dl 

FBd f ■ TmA  *-^i.-«-T-i»»d4 

PPd>  « PI  « T2di  • T3di  ♦ TMdl 

PPidl  • PI  ♦ TZldT ' 

PPZdl  ■ PI  * T22d> 


Figure  B-4  (Cont'd) . Listing  of  Shock  Propagation  Code. 
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IBIS  PASE  IS  EESI  ftUAIin  PEACH'!^ 

SSSOil  COiP^  j',UPlN'ISH£iO  rO  DD,C  ^ 


FUNCTION  UORHOKRH* 

CO'IHON  /!/  S£T*HS«QRHtTH«BrHtN«HS«Pl»B£T««|.l,LlPI»*^*SOI>l 

• i<NT 

R£*L  N.HS  ■ ■ - 

RHMl  » RH  - 1. 

R 

F . N • IRHMII  / BTH 

n « 2.  • RHHl  

T2  A 4,  • MS  * ORH 

TJ  « <N+1.»  • RHENHl  • 0R*F-»  ORH-  - • ■■• 

OORHOl  s F • ( T1  ♦ T2  - T3  » 



ENO 

SUHRilUVINE  OUTPUT 

COMMON  /!/  aiT*MS,QRH,TH,aTM,N,MS,Pl,e€TA,LltLlPl.Rt300»,OUTfDDRM  - 

* ,KNT 

OOMMUN-TJ/  P“<300l  ,PPl(3a01.PP2«300»tT2<300»,TZH30D»*T22CiOO» 

« ,T](30a)  «T4(300I 

r;»l  ttarM — — 

PRINT  loao 

PRINT  ia01<  MSf  PI  

PRINT  1102 

no  10  I ■ i*Li  •••  ■-; 

PRINT  1003,  PP(I)  .PPi  (1»  «PPZIl)»T2(I».T2ia>f  12211)  ,T3»I»  .TAtl) 

1-0- — CONTINIH 

RETURN 

1000  FORM4T<1M1a  “ THIS  PROGRAM  COHPUTES-TME  SOLUTIONS  TO  THE  EQUATION 
*SI  T3.  ’Mt  ASO  75. THESE  ARE  THE  QACH  AND  LEE  COMPLETE  SOLUTION 
*,  TH-  STRONG  SHOCK  SOLUTION  OF  BEL.  ANO  TURKOVlCH»S  SOLUTION.") 

1001  FORMAT  (IH0,“  THIS  DATA  IS  COMPUTEO  FOR  SHOCK  MACH  NUMBER  I “FT . 3/ 

tHOTttO^Pt-*“ei-2r91 

1302  FORMAT  (IMO  ,!» “PP  ( I ) ”T 24”PP  1 ( I ) -T41"i»P2  < I > “T58"T2  « D" 

♦ T »5“r21 <I) “T92“T22  <I)"T10a"T3<I)“Tt26"r4<I)~/lH0) 

100  3 FORMAT  (T3,El2.5.T20.ei2.5.I37,F12.5»T54.il2.3.T71.E12.5 

• TA«,'.12.S.T105,tl2.5-.T122itl2.5»  - 

ENO  . ..  

SUnROUTINl  PICTURE 

COMMON  71/  BET4MStORH.TH,OTM.N»MS.PliBETA.LltLlPl.RT300) .DUT.OORH 

* ,KNT 

- COMMON  72/  PP  ( 5U0  ) .PPH300)  .PP2<300)  ,T2«300».T21  1300)  ,T22  t300) 

* .T-TtOaO)  . 74(300) 

■ C?OMnON-7  -yYZOOM -7-OMI-N»p-S-fEPK>CMINfrST^2 .-PCTU.MOOtrlHt — 11 . IFtlPl — 

A IORO-R,IPLOT.IPRINT,ISTRIP.MISCtlO),>.4ST,L3LC(5),LBLCl(5). 

8 LnLK(5)  ,LnU>  ( it  3.3),LBL2(5).MOOSVM(  D.MULPLT. 

C KIN(  1,  3«  100). ZINC  1.  3.  100) .XLGTH.ZLGTH 

MODSYMCl)  ■ -1  - ‘ 

IPRINT  « 1 

VOLC  tt  )-»-18M TTtT1A-» — 

L8LC(2)  » IMA JE(0ETA,5) 

LnLC<3)  » lOH  MS  ■ - - 

L3LC(4)  > IHAGE(MS«5) 

LBLX(6)  * lOH  RADIUS — 

LBLZd)  » iOHPRiSSURE 

L-BLni-.-lTl-t-^ -KWBACN  -AWO-t 

LBLY(l,l,2)  » lOHEE 


Figure  B-4  (Cont'd) . Listing  of  Shock  Propagation  Code 


Sf  quality  FRACriQ4m 


LBl.Ytl.'iat  ■ lOMB  *HO-t 

LaLV(l«2,2i  < lOHSTHONS  SMO 

- teUV«l,2,3»  » lOHCK 

LBLy(l»3ii>  > IOHYU^KOVXCH 

tflLY<lt3te>  ■«  leHSTHOHB-SHB 

UBLYdiJiSI  « IQHCK 



C*UL  STA^T 

--BO-AO' t ■ 

XXN(l<l<t)  > XJNd,2,l>  * XIN(1«3.2)  ■ ftlX) 


XX*(a»2d>  • PPKXl 

■ f»PZ<X4 

10  CONTXHU': 

LUST  •-0 

Xf  I K>ir  .EQ.  3 ) UASX  • 1 

— e^lU-XYgPiT 

RETURN 


Figure  B-4  (Cont'd) . Listing  of  Shock  Propagation  Code, 
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